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SECTION I 
INTRODUCTION 

One purpose of thermionic-emission theory and space-charge theory 
is to relate in a sensible way the current drawn from a cathode to the 
voltage applied to the anode. Tests are made on cathodes in geometries 
kept as simple as possible to make it easier to relate measured current 
density J as a function of voltage V back to the basic characteristics of 
the cathode. The usual cathode property of interest is the cathode work 
function or the zero-field current density J q which is an exponential 
function of * and the temperature. Once the basic properties of the 
cathode are known the performance of the cathode can be predicted for 
the desired application. 

Different techniques are available for relating the zero-field 
current density J Q of the cathode to the J{V} measurements depending on the 
magnitude of J and the range of voltage used in the test. The develop¬ 
ment of theory covering three modes of operation - two which neglect the 
effect of space charge and one which includes space-charge effects - is 
reviewed below. (1) If the voltage is sufficiently large to overcome 
space-charge effects in the interelectrode space, the current becomes an 
exponential function of the square root of the voltage. This relationship 
between the current and voltage is known as the Schottky effect. (2) If 
J Q is sufficiently small, space charge effects can be neglected for 
values of J up to J Q . For this experimental condition the data do not 
show the increase of current due to the Schottky effect until the voltage 
is increased far beyond the value needed to obtain J Q . In this case the 
current is an exponential function of the voltage up to and then a 
very sharp saturation for J greater than J Q . (3) If space-charge effects 
are significant at J Q , the current departs from the exponential behavior 
at low voltages and rises to J 0 according to the space-charge theory. 
Following a review of the historical development of theory for these 
three modes of operation is an overview of the new additions to this theory. 
The new additions to the theory are then developed in detail in the body 
of the report. 


1 






AFWAL-TR-81-1019 


1. BACKGROUND 

The exponential nature of the current as a function of voltage in the 
retarding field range is discussed by Davisson (Reference 1). If space- 
charge effects are small, the current density saturates abruptly and 
remains steady at J = J . Shelton (Reference 2) gives experimental data 
which fit this theory very precisely and shows how to evaluate the basic 
cathode properties in terms of the experimental data for his particular 
example. Medicus (Reference 3) generalizes the work of Shelton and 
Davisson to develop a systematic method for analyzing retarding field data. 

Schottky (Reference 4) explains the exponential rise in current with 
the square root of the voltage for large voltages (accelerating field at 
the cathode). In the absence of space-charge effects the applied potential 
interacts with the surface fields near the cathode to lower the work 
function. The fact that the surface field is accurately represented by 
the image of the electron in the conducting surface has been established 
by the agreement of Schottky theory with experimental data (Reference 5, 
p 109). Hutchison (Reference 6, p 305) gives a review of Schottky theory. 

Langmuir (References 7 and 8) summarizes and extends previous work to 
give a theoretical treatment of the effects of space charge in the 
interelectrode space. Langmuir develops functions n and r representing 
the voltage and distance in the space between a parallei-piane cathode 
and anode and outlines a series of steps to calculate the current as a 
function of voltage for this type of diode. Van der Ziel (Reference 9) 
shows that Langmuir's n, >. functions can be extended to cover retarding- 
field and accelerating-field conditions as well as space-charge 
conditions, however, he neglects the effect of the electron image field. 

Van der Ziel mentions the transition point from retarding field to space 
charge conditions and names it the "exponential point". Although Van der 
Ziel gives a method of computing n and ' for accelerating field conditions, 
his theory still predicts that the current density is identically J for 
accelerating field conditions since he neglects the Schottky effect. 

Kleynen (Reference 10) gives several useful formulas for calculating n 
and ■' and extends the tables of r,, • values given by Langmuir. 
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Ferris (Reference 11) takes the space-charge theory developed to 

this point and presents it as a family of normalized JlV) curves which can 

be used to analyze experimental data. The experimental data are presented 

in normalized form and superimposed on the family of JM/1 curves allowing J Q 

to be determined. The normalizing constant introduced by Ferris is I, o 

which represents the maximum current density which can be drawn to the 

anode under retarding-field conditions for the condition that I tends to 

0 

infinity. Crowell (Reference 12) presents a slightly amended normalization 
technique which simplifies the family of curves presented by Ferris. 
Nottingham (Reference 13, p 65) uses I , the current at the "exponential 
point", as a normalizing constant to give a family of JiV) curves only 
slightly different from those of Ferris. Nottingham presents experimental 
data showing that the second derivative of J{V} has a sharp peak at the 
exponential point. However, Nottingham uses the second derivative only 
in conjunction with his family of JlV) curves to evaluate 0 Q . Finally 
Rittner (Reference 14) presents Langmuir’s n, f tables as interpolation 
formulas suitable for use on a digital computer and shows that Langmuir's 
theory accurately matches experimental J {VI data. 

Aside from this main development of the space-charge theory Crowell 

(Reference 15) discusses a method using the second derivative of JIV1 

in an approximate determination of J . Crowell uses a numerical study 

covering a wide range of J Q values to show that Langmuir's theory predicts 

an inflection point in the J •> V | function. This gives theoretical 

justification to the standard practice of associating the inflection point 

of experimental J{VI data with J Q (IRE standard definition of inflection 

point current (Reference 16). Crowell gives formulas for calculating 
2 2 

dJ/dV and d J/dV implying a functional relationship for J(V}. 

However Crowell does not exhibit J{V1 in functional notation. He continues 
to rely on a series of descriptive steps for calculating J<V}. 

2. OVERVIEW 

The work presented in this report extends the space-charge theory 
of parallel-plane diodes in several ways. Space-charge effects are 
included in the accelerating field range to remove the incompatibility 
between Langmuir theory and Schottky theory. A method is developed to use 
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a measurement of J g at the exponential point to calculate without 
resorting to matching a family of curves. A number of properties of 
Langmuir space-charge theory is described in much greater detail than 
has previously been done through the use of functional notation to 
express J•; VI. These details are brought out in a discussion of the fields 
and charge in the interelectrode space. As a side light, a method is 
developed to handle the edge effects in a class of parallei-plane diode 
testers. 

The effects of space-charge in the accelerating-field range are 
accounted for through the interaction of the space-charge field with the 
electron image field. The theory calculates a J-V. curve which rises from 
J = with the Child's law slope and then approaches the J•V 1 curve 
generated from Schottky theory. Since Langmuir space-charge theory 
predicts that JiV, 1 approaches J - J q with zero slope, this theory predicts 
a discontinuity in the slope of J V r at J - J Q . This discontinuity is 
pronounced in the J- V > plots for low values of J o while it occurs over 
such a narrow range of J values for large J o that it is not visible in 
the J ; V- plots. However, the theoretical discontinuity shows up as a 
sharp minimum in experimental second derivative plots (Reference 13, p 135) 
and is expected to be apparent even for large J Q . Observation of this 
discontinuity should aid in the evaluation of cathodes even in the presence 
of multiple work-function patches. An application of this extended 
theory is given to show how a common "cathode evaluation technique" lacks 
theoretical justification. The theory shows how this "evaluation 
technique" leads to erroneous values of J Q which are dependent on the 
diode spacing. (Reference 17 gives experimental data demonstrating this 
problem.) A close examination of the transition point from the retarding- 
field range into the accelerating field range (the "exponential point") 

results in a method for determining J from a measurement of the current 

o 

density at this transition point. This evaluation of J Q is accomplished 
with a single graph of J Q as a function of J g without the matching of 
experimental J'V 1 data to a family of JTV1 curves as was previously 
necessary. The theory calculates a discontinuity in the second derivative 
which shows up as a sharp peak in second derivative data measured by 
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Nottingham (Reference 13, p 135). This evaluation technique is valid for 
high-work-function cathodes and can be applied to low-work-function 
cathodes at reduced temperature. 

Previously, Langmuir space-charge theory has been used to calculate 
values of V for given values of J through a series of steps (or flow 
chart) without the use of a mathematical expression for V as a function 
of J. By exhibiting ViJ 1 as a mathematical expression, trends can be 
observed and limits and derivatives can be calculated. As an example of 
the importance of functional notation the braces : 1 are reserved in this 
report for designating "function of" (for example ulO, J VI, ln-xt). 

The use of functional notation for the ViJI relationship has led to 
the description of a number of properties of Langmuir space-charge theory 
in much greater detail, than has previously been done. The theory shows 
mathematically that there is always an inflection point for JiVi in the 

space-charge range and this inflection point coincides with a discontinuity 

2 2 

in the second derivative d J/dV when J Q is very small. The theory shows 
that J always approaches J Q with a slope dJ/dV equal to zero. The 
limiting values of several other derivatives of current, fields and 
charge density functions are evaluated at the exponential point, at the 
transition point from space-charge to accelerating field conditions, and 
at the cathode surface. 

As a side light the report gives an extension of Schottky theory 
(without space-charge) to include the effects of edges in a class of 
parallei-plane diodes. The field solution of a plane edge above an 
infinite plane is used to estimate the departure from parallei-plane 
fields near the cathode edge. For a typical example the edge effects 
introduce a 66°C error in temperature calculations based on Schottky 
theory while introducing no significant error in the determination of the 
zero field current J . 


5 
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3. THERMIONIC EMISSION 


The behavior of a free-electron gas in a three-dimensional potential 

well is sufficient to explain the basic properties of thermionic emission 

from metallic cathodes. The number density of electrons in the free- 

electron gas with momentum values lying in the interval dp dp dp is 

x y z 

given by Reference 6. 


N{p} dp dp dp. 


... dP x d V P A. 

exp {[(p x 2 +p y 2 +p 2 


)/2m - Ep]/kT} 


:d 


where 

-31 

m = 9.109 x 10 ' [kg] = mass of an electron 
k = 1.380 x IQ' 28 [J/°K] = Boltzmann's constant 
h = 6.625 x 10‘ 34 [J.sec] = Planck's constant 


T = temperature [’K] 

Ep = energy of the Fermi level [J] 


Equation 1 is the product of the number of states existing in the potential 
well within that range of momentum, times the probability that the states 
are occupied at temperature T. The wave nature of the electron is used 
to obtain the number of energy states available. The statistical theory 
of particles, in which only two particles are allowed in an energy state 
at one time, is used to determine which of the possible energy states 
are occupied at a given, nonzero temperature. The Fermi level Ep is the 
highest energy state which is filled when all lower energy states are 
filled (at 0°K). The reference level for zero energy is the bottom of 
the conduction band for this derivation. 


If the potential well has a finite depth, with wall height <|> electron 

volts above Ep, Equation 1 shows that some electrons will have sufficient 

energy to escape. For these escaping electrons, the exponential term in 

-19 

Equation 1 will be exp ie^/kTl or larger where e ( =1.602 x 10 [C]) is 

the charge of an electron. For thermionic emitters <}> is of the order of a 
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5 

few electron volts so the exponential term is of the order of 1 x 10 . 
So, for the escaping electrons, the 1 in Equation 1 can be neglected 
compared to the exponential term. Equation 1 then becomes 

N{p} dp x dpy dp z - (2/h 3 )exp{[E F - (p x 2 +p y 2 +p z 2 )/2m]/kT)dp x dp y dp z (2) 

For convenience let the normal to the surface lie in the x direction. 
Since the integral 

CD O 

; exp { - P y /2 mkT) dp u = /2miikT (3) 

for u = x or y, the number density of electrons with momentum in the 
interval dp x is given by 

N(p x ldp x = (4-rmkT/h 3 ) exp{Ep/kT} exp {-p x 2 /2mkT} dp x (4) 

If there are no fields external to the potential well, electrons 

with momentum p x can escape where 
o 

p /2m = Ep + e* (5) 

o 

Then the total current escaping from the potential well is the charge of 

the electron e times the number density of the electrons NTp^} times the 

velocity of the electrons p x /m summed over those values of momentum from 

P x to infinity: 
o 

J o = ; p (e/m) Px N <Px> dp x (6) 

x o 

The subscript "o" denotes the absence of external fields, and J is 

o 

referred to as the zero-field current density. After integration. 
Equation 6 becomes 

J o = (4n mek 2 T 2 /h 3 ) exp T-e$/kT} [A/m 2 ] (7) 

which is the Richardson-Dushman equation (Reference 6). 
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Equation 7 is more commonly expressed as 

J Q = AT 2 exp (-e<(>/kT) (8) 

where 

A = 4nmek 2 /h 3 (9) 

and A has a numerical value of approximately 120 [A/(cm 2 °K 2 )]. 

If the electrons experience an external retarding potential of height 

V $ above f> as shown in Figure 1, fewer electrons will escape so the 

current J escaping from the cathode will be less than J . The derivation 

o 

of Equation 7 applies with the use of ;■ + V s instead of Then 
Equation 7 becomes 

J = (47rmek 2 T 2 /h 3 ) exp {-e (<j>+V s )/kT} (10) 

where the zero field subscript is omitted. Equation 10 can also be 
written 

J = J Q exp i -eV s /kT] (11) 


where J Q is given by Equation 7. 

Since the Fermi level Ep is measured relative to the bottom of the 
conduction band, current will flow in an external circuit connected 
between the cathode and an anode unless the Fermi levels are at the same 
potential. So, the applied voltage (or measured voltage) of the 

a pp 

anode with respect to the cathode is just the difference in potential 
between the Fermi levels of the anode and cathode. Figure 2 shows for 
the case where the potential minimum occurs at a collecting body (or anode) 
that the applied voltage is 


( *c ~ *a } 


The derivation through Equation 11 applies to particles of charge +e, 
although electrons are conventionally assigned a negative charge. This 
conflict causes no trouble if the direction of increasing potential is 
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taken downward in Figures 1 and 2 for the development through Equation 11. 
However when measurements are involved such as V app which assume the 
negative charge for the electron, the direction of increasing potential 
must be taken upward. For this reason V s appears with a negative sign 
in Equation 12. 

For the condition that the potential minimum occurs at the anode, 
combination of Equations 11 and 12 gives 

J - J Q exp {e(V app + * c - * a )/kT} (13) 


which shows that the current is an exponential function of the applied 

voltage. This is known as the "retarding-field" or "exponential" range 

of the J-V curve. A plot of 1n{J} as a function of V gives a 

app ^ 

straight line with slope e/kT. The current J equals J Q when 



(14) 


Figure 3 shows a schematic of the I-V Curve under these conditions. 

Note that Equation 13 is actually independent of , as long as J is 
less than J , since d> is cancelled by the J term. This is shown 
through the combination of Equations 8 and 13 which yields 

J = AT 2 exp fe(V app - ^)/kT} (15) 
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SECTION II 
SPACE-CHARGE THEORY 

Figure 1 shows a schematic of the space-charge region outside an 
emitting surface. The position and magnitude of the potential minimum 
determines the amount of current which is collected at an anode spaced 
some distance away from the cathode. If the potential minimum occurs at 
the anode, the potential decreases at each point from the cathode to the 
anode. If the potential minimum occurs at the cathode, the potential 
increases at each point from the cathode to the anode. For the inter¬ 
mediate case of a potential minimum somewhere between the cathode and 
anode, retarding-field conditions occur near the cathode and accelerating- 
field conditions apply near the anode. For this intermediate case the 
electric field -dV/dx at the minimum is, of course, zero. Since, for the 
intermediate condition, the current is greatly influenced by the mutual 
repulsion of the electrons, this is known as the space-charge region of 
the J{VI curve. 

1. DERIVATION OF LANGMUIR'S DIFFERENTIAi EQUATION 

The following discussion derives the basic equations which describe 
the effect of space-charge on JiVi. Electrons with energy less than 
Ep + e; are turned back right at the metal-vacuum interface. Electrons 
with energy greater than Ep + e: enter the space outside the cathode and 
are acted on by the external potential. At a position x outside the 
cathode surface, electrons of a given energy from the momentum distribution 
of Equation 4 have lost momentum due to the barrier at the surface and 
due to the potential field outside the surface. As indicated in 
Figure 1, which shows a schematic of the space charge region outside an 
emitting surface, the net loss of energy is Ep + e$ + eV s - eVIxl at 
position x. Note that Vlxt is measured relative to the potential 
minimum. If Equation 4 is adjusted for this lost energy, it gives the 
number density of electrons as a function of momentum at position x as 
written in Equation 16 

Nip } dp = (4nmkT/h 3 ) exp iE p /kTlexp{-(p 2 /2m + E r + ef+eV -eVIxl )/kT}dp 

X X r X r j X 

(16) 
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To the anode side of the potential minimum, the lowest energy 
electrons at position x will have energy + eVlx) due to electric field 
acceleration. To the cathode side of the minimum, electrons will be 
present with momentum from zero to infinity toward the anode. In 
addition, electrons will be present with energy from 0 to eVlx) returning 
toward the cathode from the direction of the potential minimum. Then to 
find the total number density of electrons at x the limits of integration 
will be + (2meV<x') ' with the + sign taken in the anode space (to the 
anode side of the potential minimum) and the - sign taken in the cathode 
space. This integration of Equation 16 then gives 

NT V(x 11 = [(?-nkT} 3 / 2 /h 3 ] exp { -e(V s +<f>)/kTfexpieV{x)/kT) 



•[1 + Pi(eVlxl/kT) 172 }] 

(17) 


u . t 2 


where 

Piu} = (2//tT) f e dt = error function 

0 _** 

(18) 


integral 


In Equation 17 the - sign applies in the anode space and the + sign 
applies in the cathode space. 

The electric potential outside the metal is governed by Poisson's 
equation 

d 2 V(x}/dx 2 = eN/e 0 (19) 

where 

Q 

e 0 = 1/(36tt x 10 )[C/Vm] = permittivity of vacuum 

and x is measured from the position of the potential minimum. The 
combination of Equations 17 and 19 and the introduction of the normalized 
variables 

n = eVfxl/kT (20) 

and 

f = x/x-j (21) 


11 





AFWAL-TR-81-1019 


where 

X, 2 = e Q kTh 3 exp {e(V s +ri>)/ kT >/2e 2 ( 2-rrmkT) 3/2 (22) 

yields the equation 

d 2 n/df. 2 = 1/2 exp {ni[l 7 P (AT }] (23) 

This equation represents the potential distribution as a function of 
distance from the potential minimum. The normalizing constant x^ as 
more conveniently written 

= c 0 ( kT ) 3/2 / Je (M 1/? (24) 


where J is given by Equation 10. 

Multiplication of Equation 23 by dn/df; and integration with 
respect to c gives 

(dn/df. ) ^ - n 0 ^ = [1 + PiATl] exp In) - [1 + 2(nA)^ 2 ] (25) 


where h is defined as 
o 


*0 H dn/dc! r=o 


(26) 


As before, the - sign of the + applies in the anode space and the + 
sign applies in the cathode space. In order to maintain a positive 
direction for £, from cathode to anode, the square root of Equation 25 
must be taken with a negative sign in the cathode space and positive 
sign in the anode space, giving 

anode dn./d 5 = (h 2 + [1 - PiATl] exp fnJ-02(n/TT) 1/2 ]) 1/2 (27) 

space a 0 a a a 


1/2-H-1/2 


cathode d£/dn„ - -(n + [1 + Pl^nT)] exp In )~[l + 2(n /tt) lfC ~\) 

space ' c c 


(28) 


12 









AFWAL-TR-81-1019 


The subscripts on n and n are added to distinguish the anode space and 

3 C 

cathode space respectively. The reciprocal expression d;/dn c is used in 
the cathode space since the following development of the J-V curve 
requires n as a function of ; and requires ; as a function of r 

a C 

Rittner (Reference 14) has compiled a collection of accurate approximation 
formulas for the solutions to Equations 27 and 28 for A = 0. Rittner 
claims a maximum of .01% relative error within the cited range for the 
following formulas: 

0 < n c 1 0.20; ; £ = - 2n c 1/2 + .376126r, c + .0251956n c 3/2 

+ .00069163n 2 - .0009409n 5/2 (29) 

c c 

-6.156143r> -65.9123n 2 -59.4892n 3 +4.103353r, 4 

0.20 <n <3.0; ; = ---^ 

-0.2629089 -14.62018n c -49.59515r, c z -16.06387 ,^+n^ 

(30) 


3.0 < n c , r, c = -2.55389 + /~2 exp {-r» c /2> - 0.0123 exp {-r> c } 

+ (l/3/T)[(n c /Tr) 1/2 + 1] exp {-3n c /2) (31) 

0 i n < 0.05; n = 0.25 ;}- 0.047016 ; 3 + 0.012627 ; 4 -0.003820; 8 (32) 

a a a a d 3 

0.05 < n < 11.0; n, = 0.0003053192 -0.0010857509 ; + .2509922; 2 

a a a a 

-0.04648354; 3 +0.01094699 ; 4 -0.002286296 £ 5 

a a a 

+0.0003642352 E 6 -0.00004062638 ; 7 (33) 

a a 

+0.000002943653 ; 8 -0.000000123597 ; 9 

a d 

+0.0000000022744 ;, 10 

a 
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11.0 < n : n = [(' +0.5088)/!,2552] 4/3 -1.5232(f +0.5088) 2/3 

— a a a a 

+ 1.347 -0.495 (e +0.5088)" 2/3 -0.140(£ +0.5088)" 4/3 


Equation 31 comes from Langmuir (Reference 7), and Equation 34 is the 
inverse of a series given by Langmuir. Equations 29 and 32 are series 
solutions to the differential equations. Equations 28 and 27 respectively, 
taken from Kleynen (Reference 10). Rittner developed Equations 30 and 33 
by applying standard curve fitting techniques to tables of n and r values 
given by Langmuir. Langmuir generated the n, f, values by numerical 
integration of the differential equations. (Rittner gives a different 

5 

coefficient for r than the value shown here from Kleynen. The dis¬ 
cussions in ' ubsection 2 show how to verify that Kleynen's value is the 
correct one.) 

For space-charge conditions, as discussed above, the potential 
minimum occurs between the cathode and anode and dV/dx is zero at the 
minimum. For this reason one has 


dn/df ’WV° 


for space charge conditions. A plot of n(U is shown in Figure 4 along 
with the first and second derivatives. These derivatives will be useful 
below in discussions of the behavior of JiVf. Figure 4 was generated 
through the use of Equations 29 through 34 and Equations 23, 27, and 28. 
Notice that f, approaches the limiting value of -2.55389 as shown by 
Equation 31 when n c tends to «. 

2. SERIES SOLUTIONS FOR THE n, ■ FUNCTIONS 

Although Kleynen (Reference 10) does not discuss how he derived the 
series solutions for the n» ■ functions given in Equations 29 and 32, the 
following discussion gives one way of deriving the terms which Kleynen 
gives and of deriving additional terms: 

Let h 0 = 0 and let 
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then Equation 27 becomes 

du/d£ = ([1 - Ptui] exp {u 2 } - [l-2u//7"]) 1/2 /2u (37) 

From Equation 18 the derivative of P{u} is 

dP/du = (2//T ) exp {-u 2 } (38) 

so the derivative of Equation 37 becomes 

d 2 u/dr, 2 = ( (u 2 -l)[l - Plu}] exp 1 u 2 } + [1 -2u/^T])/4u 3 (39) 

Successive derivatives give 

d 3 u/d£ 3 = ([2u 4 -3u 2 +3][l-P{u}]exp{u 2 )-[3-6u//V + 2u 3 // _ jT ]) 

•(du/dc)/4u 4 (40) 

and 

d 4 u/dr 4 = (d 3 u/dr, 3 ) (d 2 u/d;; 2 )/(du/df,) 

+([4u 6 -6u 4 +12u 2 -12][l-P{u)]exp{u 2 }-[4u 5 +8u 3 -24u]//7" +12) 

•(du/ds) 2 /4u 5 (41) 

Application of L'Hospital's rule (Reference 18) to (du/dc) 2 , d 2 u/d£ 2 

3 3 4 4 

and the zero/zero terms in d u/df; and d u/d£ gives the following values 
for these derivatives at £ = 0: 

du/dc = 1/2; d 2 u/d? 2 = -1/6AT ; d 3 u/d£ 3 = 1/16 

d 4 u/dt 4 = -29/240 /T (42) 

The Taylor series (Reference 18) about £ = 0 then gives 

utf.l = a 1 e + a 2 £ 2 + a 3 £ 3 + a 4 £ 4 (43) 
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where 

a ] = 1/2; a 2 = -1/12/7 ; a 3 = 1/96; a 4 = -29/5760.7 (44) 

The coefficient a Q = u{t. = 0} is zero since the reference voltage and 
hence u is chosen to be zero at 6 = 0. 

2 

To obtain an equation of the form of Equation 32 requires n = u , 

a 

that is: 

n a = b 2 £ 2 + b^ £ 3 + b^ £ 4 + b^ £ 3 (45) 

The coefficients b n are found from a product of Equation 43 with itself. 

2 2 

b^ = a-| ; b 3 = 2a^a^-, b 4 = 2a^a 3 + a^ ; b & = 23^ + 2a^a^ (46) 

Here b Q and b-j are zero since a Q is zero. Combination of Equations 
44 and 46 gives: 

b 2 - 1/4; b 3 = -1/1277 ; b 4 = (1/2 + 1/37/48 (47) 

b 5 = -39/5760v7“ 

Equation 32 results from a substitution of the numerical values from 
Equation 47 into Equation 45. 

In order to establish Equation 29, let 

v = 77 and £ = -z (48) 

Then the derivatives of v with respect to z are very similar to those 
shown in Equations 37 through 41. In fact, the numerical values at 
z = 0 of d n v/d z n are the same as the values of d n u/d£ n with the 
variations appearing only in the sign. The values of the derivatives 
at z = 0 are shown in Equation 49: 

dv/dz = 1/2; d 2 v/dz 2 = 1/6/7 ; d 3 v/dz 3 = 1/16; 
d 4 v/dz 4 = 29/240/7 ; d 5 v/dz 5 = l/10n + 9/128 (49) 
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The Taylor series about z = 0 then gives 

v(z! = CjZ + c 2 z 2 + c 3 z 3 + c 4 z 4 + c^z 5 (50) 

whp rp 

c n = (1/n;) (d n v/dz n ) (51) 

The function vizi is then obtained by substituting £ = -z back in 
Equation 50: 

vtr.l = g^ + g 2 c 2 + g 3 £ 3 + g 4 £ 4 + g & £ 5 (52) 

where 

9 1 = -1/2; g 2 = 1/12/T; g 3 = -1/96; g 4 = 29/5760/7" (53) 

g 6 = -(1 /IOtt + 9/128) /120 

A comparison of Equations 52 and 53 with Equations 43 and 44 show that 

ufrl = - v {r.} (54) 

where the equality is taken in the sense of the Taylor series representation. 
Consequently the square of Equation 54 shows that the Taylor series repre¬ 
sentation for n a {f,> is exactly the same as the Taylor series representation 
of n.lK) at r. = 0. This shows that the function n, taken as an extension 
of the function r> is continuous at £ = 0 and has continuous derivatives. 

To obtain an equation of the form of Equation 29 requires zTv> rather 
than vizi, that is 

z = djv + d 2 v 2 + d 3 v 3 + d 4 v 4 + d 5 v 5 (55) 

The derivatives of the inverse function z{v} can be found from formulas 
relating the derivatives of the inverse function to the derivatives of 
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the original function, some of which are given by Pierce (Reference 19, 
p 100). The formulas necessary in this case are 


dz/dv = 1/v (56) 

d 2 z/dv 2 = -v/(v) 3 (57) 

d 3 z/dv 3 = [3(V) 2 - 'v'v]/(v) 5 (58) 

d 4 z/dv 4 = [10 v v V - 15 (v) 3 - (v) 2 v^ v ]/(v) 3 (59) 

d 3 z/dv 3 = 10[( v ) 2 (v) 2 + 15 v v lv (v) 2 - 105 (v) 2 v v 

+ 105 (V ) 4 - v V (0) 3 ]/(v) 9 (60) 

where, for example, d 3 v/dz 3 has been shortened to'v, and d 4 v/dz 4 
has been shortened to v 1V . Using the formula 

d n = ( 1/n! )(d n z/dv n ) (61) 

with Equations 49 and 56 through 60 gives: 

d| = 2 ; d 2 = - 2/3.T" ; d 3 = - (1 - 8/3 tt)/6; 


d 4 - -2 (5/3 tt - 21/40)/9.> ; d 5 -- (1/2 -72/5^ + 1 120/27 t, 2 )/120 (62) 

Substitution of Equation 48 into Equation 55 gives 


r _ A 1/2 
* ~ - d l n c 


■ .3/2.2. 5/2 

' d 2\ ' 3 ’c ’ d 4 n c ' d 5 n c 


(63) 


Finally, the use of Equation 62 in Equation 63 gives the series shown in 
Equation 29. 


For the case where n = 0, Taylor series do not exist for n If.) or 

O d 

£, ■ d 1 • This is seen on examination of the third derivative of Equation 27 
d n / d f 3 = (1/2)[(1 - P {» 11 i) exp ini - T / v^rTn ] 

o add 

•[n 2 + (1 - p i,e t) expin 1 -(1-2,^)] 1/2 (64) 

U d d d 

The term 1/ > r rT increases without bound as n approaches zero, showing 

d d J 

that the third derivative of n if.l is unbounded at f, = 0. Equations 23 
and 27 show that at i = 0 

dn /df, = 6 and d 2 n,/d6 2 = 1/2 (65) 

d U d 
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The calculation of the derivatives as in Equation 69 and their 
evaluation at u = > =0 gives 


v 1 = 0; 


V 2 = 17 V 


» 3 - 0; v„ ■ -1/4^; 


(71) 


»5 ' 4/1 5- m 3 0 , 


(2/n 4 -3/.'.*)/24; » ? " 2(4/^ - 15/^)/10b/7 


The combination of Equations 67 and 6B and omission of the zero terms gives 


v 2 n a + v 4 n a + v 5 


5/2 3 7/2 

1 + V/- n + V -7 n 

3 o d id 


(72) 


Figure 5 shows a plot of n. versus c generated from Equation 72 with 

a 

n Q = 0.5. The series solution with up to the seventh power of u from 
Equation 72 is shown along with the series solution truncated at the 
sixth power of u. The full series solution is useful up to c equal to 
about 0.6. Since the coefficients involve reciprocal powers of n in 
Equation 71, one might expect the series to behave better at larger 
values of n Q . Although the series is able to handle larger values of 
n. as 6 becomes large, the range of accuracy for c is still 1. Note 
that n is a steeply rising function of r, for large values of n . The 

3 O 

function n fr] calculated by numerical integration is shown in Figure 5 
for comparison with the series solution. The numerical-integration solution 
procedure will be discussed in the next section. 

3. NUMERICAL SOLUTIONS FOR THE n, £ FUNCTIONS 

If f, = 0, integration of Equation 28 for tin ) cannot be started 
at n = f. = 0 since d f /dn c is unbounded at this point. Most numerical 
solution techniques for differential equations, such as the Runge-Kutta 
methods, can be used, however, if Equation 29 is used to start the solution 
procedure. The function njL I also has starting problems with Runge- 

a 

Kutta techniques. 

The problem with n,(U is that dn /d£, is a specified function of n 

3 3 3 

starting at n =0. Any solution technique such as the Runge-Kutta types 
which extrapolate values of n 3 based on dn /dc will always give n,td = 0 

3 3 3 

for all c if it begins with n a (0} = 0. This problem is solved through the 

3 
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use of a Taylor method of second order or higher. A method which has 
been very useful is the mixed Runge-Kutta-Taylor method (Reference 20, 

3 

p. 275). This method yields convergence of order h , where h is the step 
size, and uses the first and second derivatives of n from Equations 23 

d 

and 27. This solution method is especially convenient since the second 

derivative from Equation 23 has a simple form. The mixed Runge-Kutta- 

Taylor method for n generates a sequence of numbers n , for a fixed 
3 n 

value of h, by the recurrence relation 

o Q = 0 (73) 

n n+l = ”n + hG 3 tn n ;h} 

wi th 

G,u ; hi = fir, } + (h/2) f{n + hf{n }/3) (74) 

o 3 3 da 

where 

f = d n /dr (75) 

a 

and f and f = df/dr, are given by Equations 23 and 27. 

The accuracy of the Runge-Kutta-Taylor method solving for n Jr,} is 

a 

limited less by the step size h than it is by the accuracy of the values 
calculated for the error function Piul. The most satisfactory representation 
for P{u1 used in this work comes from Abramowitz (Reference 21, page 299, 
Equation 7.1.26). This representation is 

P{ul = 1 - (a^t + + a^t^ + a^t^ + a^t^) exp {-u^l (76) 

where 

t = 1/(1 + pu) (77) 

and 

p = .3275911 a 1 = .254829592 

a 2 = -.284496736 a 3 = 1.421413741 (78) 

a 4 = -1.453152027 a 5 = 1.061405429 
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Since the machine available for calculations in this work had only a six 
digit capability, the numerical values in Equation 78 were rounded to 


p = .327591 

a 3 = 1.42141 


a 1 .25483 a 2 = -.2845 

a 4 = -1.45315 a 5 = 1.06141 


(79) 


It is important in this round-off procedure to have the terms through 
a^ to add to exactly 1. For example the rounding of a 2 to -.284497 instead 
of the value in Equation 79 results in greatly increased error in the n 

d 

values near e, = 0. The use of Equation 76 for P{u> along with the 
numerical values in Equation 79 in the Runge-Kutta-Taylor solution 
procedure generates accurately the values of n {e,> tabulated by Kleynen. 
With this representation of P{ul, step size h less than h = .01 does not 
increase the accuracy of the solution. Step size h = .01 generates the 
Kleynen (Reference 10) tables for n {r,} accurate to the fourth place 
for f, ^ 5, and h = 1.0 gives good accuracy above K = 100. A BASIC 
computer language program using Equations 73 through 79 for solving 
Equation 27 is shown in Figure 6. A sample print-out from the program 
for h = .01 is shown in Figure 7. 

For the case where f, is greater than zero a first examination seems 
to indicate that a Runge-Kutta method would be better for solving for 
n a {} than the Runge-Kutta-Taylor method of Equation 74. Since the third 
derivative. Equation 64, of o l£} is unbounded at f, = 0, one might 

a 

expect problems with the calculation of » from Equations 73 and 74 which 

a 

use values of the second derivative of n . Note that 

d 

fin 1 = d Z n Jdr 2 (80) 

a a 


This concern is based on the Taylor Theorem. When a function is 
generated from its first and second derivatives, the remainder term is 
related to the third derivative. In addition, when n is greater than 
zero, the Runge-Kutta methods do not generate the spurious solution 
n. = 0 as discussed earlier. 

d 
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However, in practice, the Runge-Kutta-Taylor method yields much 
quicker convergence than the simplified Runge-Kutta method which uses 
the same number of function evaluations per step. This superior conver¬ 
gence is obtained over a wide range of values of n . The simplified 
Runge-Kutta method is described by 


n 


o 


= 0 


wi th 


Vl = n n + hK 2 { V h} 


K_in ;h> = (l/2)[f{n_ > + fln a + hfln M] 

c a a a a 


(81) 

(82) 


where 

f = dn /dr, (83) 

a 

and f is given by Equation 27. The convergence of the Runge-Kutta- 

Taylor method compared to the convergence of the simplified Runge-Kutta 

method is shown in Table 1 for the example h Q = .003 at 5 = 5.0. Later 

work will show that these conditions correspond approximately to a 

cathode with I 30 mA/crn^, T = 650°C, I = 1.0021 o , and cathode-anode 

spacing 0 = .05 mm. Table 1 also compares the convergence of these two 

methods for the example n = .5 at £ = 1000. Later work will show that 

o c 2 

these conditions correspond approximately to a cathode with I = 5A/cm , 

T = 900°C, 1=1.11, and 0 = 1.0 mm. 

o 

The superior convergence of the mixed Runge-Kutta-Taylor method 

3 

results from its third order convergence, order of h , compared to the 

2 

second order convergence, order of h , for the simplified Runge-Kutta 
method. For small values of the simplified Runge-Kutta method 
generates tne solution n = n r which has nearly the error of the spurious 
solution n = 0 generated for the case h =0. On the other hand, the 

3 O 

third derivative of n is large over a very small interval of values of £ 
for small h . For this reason, the total change in the second derivative 
is not large near £ = 0 and the Runge-Kutta-Taylor method accurately 
calculates n a • For large values of h 0 the value of n a is dominated by the 
first derivative so that large variations in the second derivative 
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TABLE 1 

COMPARISON OF THE CONVERGENCE OF THE MIXED RUNGE-KUTTA-TAYLOR, R-K-T, 


METHOD WITH THE 

SIMPLIFIED RUNGE- 

KUTTA, R-K, METHOD 

n Q = 1 x 10 -3 £ = 5 



step size 

n from R-K-T 
a 

n from R-K 
a 

.005 

3.61135 

3.61046 

.01 

3.61139 

3.60858 

.05 

3.61163 

3.58151 

.1 

3.61215 

3.53811 

.5 

3.63247 

3.13272 

1.0 

3.69794 

2.62286 

5.0 

5.79598 

.300455 

n Q = 0.5 £ = 1000. 




step size 

n from R-K-T 
a 

from R-K 
a 

.5 

7287.76 

7286.99 

1 

7288.19 

7286.38 

5 

7291.94 

7276.32 

10 

7297.70 

7259.18 

100 

7494.38 

6838.45 

200 

7804.92 

6328.12 
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near 5 = 0 do not cause problems with the accuracy. Consequently, the 
mixed Runge-Kutta-Taylor method accurately calculates n a also for large 
h . Since the mixed Runge-Kutta-Taylor method converges faster for 
small n Q and at least as fast for large h Q it is used for the calculations 
reported later in this work. The criterion for the choice of h, used for 
n Q >0 is the same as the criterion for h Q - 0 discussed above. 

4. THE CURRENT - VOLTAGE RELATIONSHIP FOR J <_ J Q 

If J < J , Equations 27 and 28 give the voltage according to the 
procedure outlined by Langmuir (References 7 and 8 ). If the space- 
charge potential minimum occurs at the anode, the voltage can be 
calculated from Equation 13. Otherwise, the potential minimum occurs 
somewhere between the cathode and the anode, and the value of the 
potential minimum relative to the metal surface is V $ from Equation 11: 

V $ = (kT/e)ln (84) 

The distance from the potential minimum to the cathode is found from 
Equation 28 since at the cathode n c = (e/kT)V s : 

£ = £ in ) = £ r Un{J /J}} (85) 

La L U 

cathode 

Recall that £ c and n c are measured relative to the potential minimum and 
x = £ x^ from Equation 21. The distance from the potential minimum to 
the anode is then the anode - cathode spacing D minus the distance from 
the potential minimum to the cathode given by Equation 85. In normalized 
terms, then, the locations of the anode is £, where 

Q 

£ = 0/T/ x„ + £ r { 1 n{J /J}} ( 86 ) 

a CL O 

Here the normalizing constant x^, based on Equation 24, given by 

x 2 2 = e Q (kT ) 3/2 /e( 2 *m) 1/Z = x^J (87) 
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allows the current density J to appear explicitly in Equation 86. Then 
from Equation 27 the potential V a of the anode relative to the potential 
minimum is 

V a = ( kT /e)n a {DAT /x 2 + c c tlntJ Q /J}l (88) 

Figure 8 shows a schematic of the potential distribution between 
the anode and cathode under space-charge-1imited conditions. This 
diagram shows that the anode to cathode voltage V is 

V = V - V. (89) 

d b 


and the applied voltage which is offset by the Fermi levels is 


V 

app 


V + * a 


V 


C 


(90) 


(A similar situation is discussed above in the derivation of Equation 12.) 
Then from Equations 84, 88, and 89 

V = (kT/e)[ln{J/J Q } + n a (D /T /x 2 + f, c ( ln{ J Q /J}} ] (91) 

where h Q = 0 in Equations 27 and 28 and D is the cathode-anode spacing. 

Equation 91 is instructive as a summary of the solution procedure 
outlined by Langmuir. Each term in Equation 91 represents a part of the 
physical behavior of the emission. The term ln{J/J Q } is negative and 
represents a voltage offset due to the existence of the potential 
minimum V . The term r {1n{J /JH is negative and represents the 
distance from the potential minimum to the cathode. The term Dv'T” /x^, 
serving as the argument for the function ri , is responsible for the 
basic 3/2 power behavior of the current as a function of voltage. 
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5. BEHAVIOR OF J{V} FOR SPACE-CHARGE-LIMITED CONDITIONS 

Three characteristics of the J{V} curve are of special interest for 

space-charge-limited conditions. Each of these three characteristics, 

the transition from retarding-field to space-charge-limited operation, 

the inflection point, and the transition from space-charge-limited 

operation to the zero-field condition J = J , can be used in an 

experimental determination of the zero field current density J which is 

o 

characteristic of the cathode under test. The derivative of Equation 91 


dV/dJ = (kT/eJ)[1 + (D/3~/2x 2 - d£/dn c 


) dn/df; 
a 



(92) 


and the second derivative 


d 2 V/dJ 2 = (kT/eJ 2 ) 


-1 - (dn /dr, 

a 


) [-dc/dn. 


a 


+ ST D/2x 9 ] 

I C. 


? ? 
+ (d n /df, 

a 


f , 1 

a 


) [-df,/dn f 


+ ST D /2 x 2 Y 


(93) 


+ (dn / df, 

d 

) [ - (dc/dn 

r 1 L 

) 3 (d 2 n /dr 2 

1 v. 

) + D/J“/4x 2 ] 


r 

a 

n c 

n c 


show more explicitly the behavior of the current-voltage curve for 
space-charge-limited conditions. Here n' and £ 1 are defined as 

C 3 


and 


= ln{J Q /J} 

(94) 

= ST D/x 2 + L c (ln{J 0 /J}} 

(95) 


2 2 

The derivatives dJ/dV and d J/dV are obtained from Equations 92 and 93 

through the use of Equations 56 and 57. Equation 57 is also used to 

2 2 

express the derivative d f./dn in Equation 93. 
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Figure 9 shows J{V}/J Q , (dJ/dV)/J Q and (d^J/dV^)/J Q for an example 
cathode having J Q = 10 mA/cm^, T = 700°C, and D - 0.2 mm. The curves are 
normalized by the factor J Q for ease of plotting. The curves in Figure 9 
are shown dashed in the retarding field range and solid in the space-charge 
range. The retarding-field range is described by Equation 15. 

At the transition from the retarding-field range to the space-charge 

range. Figure 9 indicates that J{V} is a smooth curve and that dJ/dV is 

2 2 

continuous. The second derivative d J/dV , however, is discontinuous at 
this transition point. The second derivative is much larger on the 
retarding-field side than it is on the space-charge side. (In the 
generation of Figure 9, Equations 92 and 93 will generate very curious, 
though spurious, values when applied in the retarding-field range.) 

Toward increased voltage in the space-charge range, notice that 
2 2 

d J/dV passes through zero, indicating an inflection point in the JfV} 

curve. Crowell discusses the existence of this inflection point through 

plots of the results of many numerical calculations. However he does not 

relate the behavior of the current-voltage curve back to the equation for 

J{V}. Crowell also does not point out that for practical, high-current- 

density cathodes this inflection point occurs virtually at J = J Q . Note 

that the example in Figure 9 is for a cathode with very low current 
2 

density, 10 mA/cm , and a very close diode spacing of 0.2 mm. A detailed 
discussion of this inflection point follows below in Subsection b. 

Perhaps more important than the existence of the inflection point is 

the fact that J{V} tends to J Q with zero slope. Work discussed below in 

Subsection a will show that this behavior of J{V) near J q is characteristic 

of high-current-density cathodes as well as the low-current-density 

cathode of Figure 9. The work below and the discussions in Section IV 

on the behavior of J{V} in the accelerating-field range predict 

2 2 

discontinuities in dJ/dV and d J/dV at J - J Q . This behavior should 
show up as rapid variations in these derivatives for cases where experi¬ 
mental conditions (such as edge effects or emission patches) obscure 
details of the J{V> curve. These rapid variations in the derivative 
should serve as a better indicator of the value of J Q than the inflection 
point measurement discussed by Crowell (Reference 12). 
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a. Transition to J = J Q 

At the point where J is just equal to J Q from the space-charge 
range, the space-charge minimum occurs right at the cathode. In this case 
n’ from Equation 94, which represents the voltage offset due to the 
potential minimum, is zero. Then in Equation 92 since dn c /dc is zero at 
this point, the term dt/dn c is unbounded. Then dV/dJ is unbounded which 
implies that 

dJ/dV = 0 at J = J Q (96) 

This is true for any value of J Q . This behavior appears clearly in an 

example treated by van der Ziel (Reference 9). (See van der Ziel's 

2 

Figure 1 with I $ = 0.01 A/cm .) 

For many practical problems this behavior of the space-charge 
equation only applies over values of J extremely close to J Q and is not 
readily visible in the current-voltage plots. For example, the zero slope 
of J as a function of V at J = J Q is not visible in Figure 24 which 
treats a high-current-density cathode. However, later work in Section IV 
shows that dJ/dV is greater than zero at J = J Q from the accelerating- 
field range. This discontinuity in the first derivative and the accompany¬ 
ing discontinuity in the second derivative should show up in plots of 
experimentally taken derivatives even when the transition at J = J Q is 
masked in the J{V} plot itself. 

b. The Inflection Point 

The inflection-point current is defined as an approximate 
measure of J Q in the 1950 I.R.E. Standards on Electron Tubes (Reference 16). 
This standard is based on the observation that J{V} has an increasing 
slope for small J and a decreasing slope for large J. The inflection 
point is then the transition point from increasing to decreasing slope of 
J{V>, the point where d J/dV is zero. Crowell (Reference 12) was the 
first to point out that Langmuir space-charge theory predicts an inflection 
point for some J less than J Q . Crowell used a numerical study for many 
values of J Q to support the use of the inflection point as an approximate 
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measure of J Q . This section goes beyond the work of Crowell by showing 
on a theoretical basis how Langmuir space-charge theory always produces 
an inflection point at some J less than J Q . 

Equation 93 along with Equations 92 and 57 give a fairly 

2 2 

complicated expression for d J/dV . However, in the limit as J approaches 
J , the term df./dn is unbounded. Use of this fact allows a simple 
expression for d^J/dv at J = J Q 


2 2 , 

d J/dV 

| J 


= -(J/2)(e/kT) 2 (dr, a /d; | ) 


(97) 


At J = J fl , f, is zero so from Equation 95 

f' - /J” D/x 0 

a c 


(98) 


Then Equation 97 shows that the second derivative of J{V) is always 

2 2 

negative and finite at J . Although Figure 9 indicates that d J/dV has 
a large negative value for that example. Equation 97 shows that the 
second derivative is finite. 

For small J, in the retarding-field range. Equation 13 shows 
that J{V} and its derivatives are 

J = J Q exp {e (V app + <fr c - <}> a )/kT} (99) 

dJ/dV = (eJ o /kT) exp le (V app + ; c - D a )/kT} (100) 

and 

d 2 J/dV 2 = (e/kT) 2 J o exp (e (V app + ^ - * a )/kT) (101) 

Equation 101 shows that the second derivative is always positive for 

2 2 

small J. Comparison of Equations 97 and 101 shows that d J/dV changes 

sign for some value of J between J = 0 and J = J . In other words, 

o 

there is an inflection point of JfV} at this value of J. 

This inflection point has the usual properties except for the 
case of very small values of v^J - O/x^ discussed in Subsection c. In this 
case a discontinuity in the second derivative coincides with the inflection 
point. The second derivative then passes discontinuously from positive 
values to negative values. 
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c. Transition from Retarding-field to Space-Charge-Limited 

Conditions 

This section explores the detailed behavior of Langmuir space- 
charge theory at the transition from the retarding-field range into the 
space-charge-limited range. Van der Ziel (Reference 9) mentions this 
transition point and gives it the name "exponential point". However 
van der Ziel does not discuss the behavior of the derivatives of J{V} 
at the exponential point, probably since he did not have Equation 91 or 
an equivalent. The examination below of J{V] and its derivatives show 
that the J{V> curve is smooth at the transition from the retarding field 
region into the space charge region ("exponential point"). In addition, 
the derivative dJ/dV is continuous but not smooth, and the second derivative 
is discontinuous at the exponential point. 

For the condition that the space-charge just begins to reduce the 
current from the value it would have in the retarding-field range, the 
space charge minimum occurs right at the anode. In this case the dis¬ 
tance from the potential minimum to the anode is zero so the value of £ 
at the anode, designated r' above, is zero. From the definition of £' 
in Equation 95, it follows that 


c c nn(J o /J e n - - ^ d/x 2 


( 102 ) 


where J is the current density at the exponential point. Then since n , 

c d 

evaluated at r = 0, is zero, Equation 91 reduces to 


V e = (kT/e) In {J e /J Q } 


(103) 


where V g is the voltage at the exponential point. Since h , which is 


dn/d£ evaluated at £ = 0, is zero, Equation 92 reduces to 


dV/dJ! 


J = J 


= kT/e J, 


(104) 


and Equation 93 reduces to 

d 2 V/dJ 2 i = (kT/e J ' 2 )[ -1 + 1/2 (-d£/dn. + /J~D/2x 9 ) 2 ] (105) 

J = J 6 C n' 

1 e c 

2 2 

In Equation 105 use has been made of Equation 23 to evaluate d n /d£ at 

3 

^ = °- 
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If Equations 103 through 105 are written in terms of J{Vj instead 
of V{J} they become 

J e = J Q exp { eV /k T} ( 106 ) 

dJ/dV - ej /kT (107) 

• J = J e " 
and 

d 2 J/dV 2 = (e/kT) 2 J [1 - 1/2 (-d ./&. , + .TO/2x ? ) 2 ] (108) 

'j = J e ^ 

e 1 c 

Comparison of Equations 99 and 100 with Equations 106 and 107 shows 
that J{V} and dJ/dV are both continuous at the exponential point and 
consequently J{V} is smooth. Note that V and V - :• are the 

d pp C a 

same from Equation 90. Comparison of Equations 101 and 108 shows that 
2 2 

d J/dV is always greater on the retarding field side of the exponential 
point than it is on the space-charge side. From Equation 28 the term 
d.-'./dn in Equation 108 is always negative. (This is also seen in 
Figure 4.) Then the term in brackets in Equation 108 is always less 
than the value calculated from Equation 101. Then the second derivative 
is discontinuous so that the first derivative is not smooth at the 
exponential point. Figure 10 shows a detailed sketch of dJ/dV at the 
exponential point for the example of Figure 9. 

For many cathodes where J o is relatively large, ■' c l ln{ J Q /J 1 1 
is very nearly its minimum value of -2.5539 before it is equal in 
magnitude to D»"J /x^. In this case J g at the exponential point is 
much less than J Q and, from Equation 102, J g is given by 

J e * (2.5539 x 2 /D) 2 (109) 

In this case of large J Q then, at the exponential point becomes 
nearly independent of the zero field current J Q . In this case dt/dn 
at is very small, so the second derivative is approximately 

d 2 J/dV 2 : -- (e/kT) 2 J [1 - (l/2)(?.5/2) 2 ] (110) 

1 1 - 1 “ 
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which is approximately a factor of 5 jump in the second derivative at the 
exponential point. For the example of Figure 9 the second derivative 
jumps by a factor of 5.7. Note in Figure 9 that the second derivative 
is off scale to the retarding field side of the transition point. 

For cathodes with very small J Q the term AT D/x^ in Equation 102 
will be very small leading to a very small K '• In this case -dt/dr 

c 2 2 ^ 

would be quite large in Equation 105 rendering a 3/6 v negative. If 

2 2 

d J/dV is negative to the space-charge-1imited side of the exponential 
point, there is no additional inflection point in the J{V} curve and the 

slope of J{V} always decreases in the space-charge-1imited range. Since 

2 2 

d J/dV is positive in the retarding field range, this discontinuity to 

2 2 

negative d J/dV can be quite pronounced for low-current-density cathodes. 


d. Experimental Use of the Exponential Point 

The discontinuity in the second derivative at the exponential 

point should lead to a sharp peak in the second derivative of experimental 

J{V) curves. The current J and the applied voltage 1/ „ can then be 

e r app 

measured at this peak in the second derivative. Then, if J Q is not too 

large, the theory developed in Subsection c and expanded below calculates 

the zero field current J and the work function difference 4> - 4> - 

o a c 

A suggested limit on J Q for accurate calculations is that be greater 
than approximately J Q /50. This condition on J Q is satisfied for large 
tj> cathodes or by sufficiently low temperature T. 


Based on one measurement of current J and voltage V , the zero- 

e e 

field current J Q of the cathode and the difference in cathode and anode 

work function can be calculated without reference to the Richardson 

Equation 8. Equation 8 can then be used to calculate the cathode work 

function from the zero-field current density. From Equation 90 the 

measured voltage V, „ is related to the voltage V„ by 
app 3 e J 


V = V 
app e 



( 111 ) 


From Equation 102 the current J g is related to J Q by 

J o = J e exp{n c {,/J e d/x 2 }1 


( 112 ) 
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and Equation 103 relates V to J and J . Then Equations 103 and 112 can 

be used to generate plots of J Q and V g as functions of J g . For a measured 

value of J & and with measured values of T and D, the zero field current 

density can be calculated from Equation 112. Then Equation 103 can be 

used to calculate V . From the calculated value of V and the measured 
e e 

value of V , Equation 111 then produces the difference between the 
app r 

cathode and anode work functions U - ). 

C 3 

The problem with the use of Equation 112 is that it requires 
n c ft) in the cathode space. A method for calculating this function could 
be derived in a straightforward manner. However, since the inverse 
function } can be calculated from Equations 29 through 31, it seems 
easier at this point to restate the problem, so that this inverse function 
can be used. The use of this inverse function is possible through the 
construction of an auxiliary function f{J} where based on Equation 102 

f{J} = £ c (ln{J 0 /J}} + 0/T /x 2 (113) 

Then finding the solution J to Equation 112 for a given J Q is equivalent 
to finding the solution to Equation 113 such that 

f{J e ) = 0 (114) 

This allows a value of J g to be generated for each J Q over a range of 
values of J Q . This can, of course, still be plotted as J Q {} for use 
in interpreting experimental data. 

A plot generated from the BASIC program discussed below is shown 
in Figure 11 for two examples. In one example the temperature is 700°C 
and the spacing is 0.2 mu while in the other example the temperature is 
900°C and the spacing is 1.0 mm. Actually, the values of temperature and 
spacing are not unique to each curve in Figure 11. Rather, the family of 
curves suggested in Figure 11 depends on the parameter t where 

o z A 1/2 

t = D/x 2 = (D/T 3 / 4 )(e/2irm /c Q k 3/2 ) (115) 

as seen in Equation 113. Figure 11 also indicates the asymptotic values 
where J g ceases to increase for increasing values of J . The asymptotic 
value is, of course, dependent on cathode temperature and diode spacing 
and is given by Equation 109. 
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Equation 114 can be solved for J g by an iterative technique 
known as Newton’s method (Reference 20, p 277). Based on an initial 
guess J-j at the solution, this method generates a sequence of approximations 
by the recurrence relation 

J n + 1 = J n - f{J n )/f{J n } n = 1,2,3. (116) 

The iteration is terminated when J , is sufficiently close to J and 

J . . is taken as the value of J . Since J is always less than 
n + 1 ^ e e J 

(2.554 x^/D; from Equation 109, this is a convenient starting value for 

J,. However, for small values of J this formula gives a value of J 
1 o e 

which is greater than J Q . In this case a smaller starting value for J-j 
is necessary, for example J-j = 0.5 J . In Equation 116 the notation 
f' {J^l represents df/dJ evaluated at J , that is 

f'{J n ) = (-l/J)(df./dn c ) + D/2 x 2 , / 3T (117) 

where df/dn is calculated from Equation 28 with n evaluated at 

n c = ln{J Q /J n > (118) 

A BASIC program which solves Equation 114 and generates and V g for a 
range of values of J Q is shown in Figure 12. 

Another way of representing the behavior of J{V} at the 

exponential point is shown in Figure 13. This figure shows both 

and J e plotted as functions of V . The calculations were made by the 

BASIC program of Figure 12. Figure 13 is especially convenient for 

interpreting measured data since it gives both J Q and V 0 based on the 

measured value of J . From V and the value of V „ at the transition 
e e app 

point. Equation 111 then gives the difference in work function ( 4 > - <j> ). 
Figure 13 shows the asymptotic value of J g for the example T = 700°C 
and D = 0.2 nun. Figure 13 also shows the asymptotic line for J Q (V e } 
which passes through the asymptotic value of J g at V p = 0 with a slope 
of kT/e. The existence of an asymptotic line for constant is seen 
from Equation 103. For experimental interpretation families of curves of 
J Q or J g as functions of eV/kT can be generated just as the sample curve 











AFWAL-TR-81-1019 


was generated for Figure 13. With the V axis normalized to eV/kT, the 
family of curves is generated with the single parameter t = D/x^ just as 
in the generation of the family of curves like those in Figure 11. 

Although the representation is not as intuitive as the families 

of curves suggested in Figures 11 and 13. All of this information can be 

represented in a single plot of .1^{V e } and a single plot of J Q {} through 

proper normalization of the terms. The proper normalization is 
2 2 

t J and x J as functions of eV /kT. Graphs of these two normalized 
e o e 

functions are shown in Figure 14. 

The derivation for these normalized functions follows from the 
combination of Equations 102 and 103 to give 

J e = (x 2 /D ) 2 C c 2 {eV e /kT) (119) 

The combination of Equations 106 and 119 then gives 

J o = ( x 2 /D)2 f ' c 2{eV e /kT} exp { - eV e /kT} (120) 

The normalized functions are then from Equation 115. 

t 2 J 6 = C c 2 { eV e /k T ] (121) 

and 

t 2 J q = T c 2 (e V e /kT) exp {-eV e /kTl (122) 

The function f, is conveniently calculated from Equations 29 through 31. 

c 2 

Figure 14 shows the asymptotic value of 6.5224 for t J and the 

2 e 

asymptotic slope of 4 for t J as a function of eV/kT near V = 0. 

The limiting value of t is just the square of the limiting value of 

S c from Equation 31, and the limiting slope is apparent from the square 
of Equation 29. 

Equation 112 gives J Q and consequently the cathode work function 
(from Equation 8) based on the measurement of current, cathode temperature, 
and diode spacing. It does not require a knowledge of the anode work 
function or the value of the applied voltage. Equation 112 in principle 
allows the determination of J Q for any value of J g . However, for large J Q 
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the value of /J~~ D/x^ has nearly the limiting value of -2.5539. In this 
case n c is such a steeply rising function that the method is not useful. 

In practice the problem with high values of J Q can be overcome by lowering 
the temperature, since J Q is an exponential function of the temperature. 

6 . COMPARISON OF J{V} WITH CHILD'S LAW 

For many practical cathode tests the temperature is around 1300°K, 

the cathode-anode spacing is of the order of 1 mm and the current density 

2 

J Q is of the order of 1 A/cm requiring 100V or more to obtain J Q from 
the cathode. Then in Equation 91 the term kT/e is of the order of .1 volt, 
and the voltage offset due to the potential minimum, (kT/e)1n{J/J q } , is 
about 0.2 volts for J = J Q /10 and decreases as J approaches J . This 
term, then, has very little effect on an applied voltage in the 100 volt 
range. An examination of other terms in Equation 91 shows that D/J - /x^ 
is in the range of a few hundred. Since £ c has a maximum magnitude of 
2.5 and decrease for J approaching J , this term has only a small effect 
on the argument of n a in Equation 91. Then Equation 91 is accurately 
approximated by 

V app ' (kT/e) V DAF /x 2 } (123) 

The function n (U for large f, can be approximated using Equation 27. 

a 

For large r., n fl is also large as seen for example in Figure 4. For large 
n g Equation 76 shows that the function P{u> is accurately approximated by 

P{u} 5 1 - exp {-u 2 } (124) 

Then for large £, Equation 27 is accurately approximated by 

dn /df, = (2/ir 1/4 )n 1/4 (125) 

a a 
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Figure 15 gives a graphical comparison of the approximate expression 
from Equation 125 compared to the exact value of the derivative calculated 
from Equation 27. This graph shows that the magnitude of the error in 
using Equation 125 decreases slowly as f increases. However the error 
expressed as a percentage of the value of dn /dt becomes small very 

a 

rapidly for c greater than 50. In generating Figure 15, Equations 32 
through 34 were used to calculate n U). Tnen Equation 27 and 125 were 

a 

used to calculate the derivative and the approximate derivative as 
functions of f;. Integration of Equation 125 shows that 

n a (U = (3/4)(3M 1/3 C 4/3 (126) 

Here the integration constant is zero since n (0) is zero. Equation 126 

a 

also follows from Equation 34, though less intuitively. The number 
1.2552 in Equation 34 is just [(3/4)(3/tt) I / 3 ] 3 / 4 rounded to 5 digits. 

The combination of Equations 123 and 126 gives 

V app = [(9/4)(0 2 /e o )(m/2e) 1/2 j' 2/3 J 2/3 (127) 

or, which is equivalent 

J = (4/9)e Q (2e/m) 1/2 V 3/2 /D 2 (128) 

Equation 128 is known as Child's law (Reference 7). Review of the above 
procedure then shows that the term 2/n /n in Equation 27 is responsible 

d 

for the basic three halves behavior of J{V}. 
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SECTION III 
SCHOTTKY THEORY 

The space-charge-theory developed to this point in Section II suggests 

that there will be a zero-field current density J o for some voltage V Q and 

that there will be no increase in J above J for V > V . However, 

0 0 

Schottky (Reference 4) suggests that a linear applied accelerating field 

across a diode will interact with the intrinsic field near the cathode 

surface to lower the cathode work function. Consequently, the current 

density will increase above J for accelerating fields at the cathode 

o 

surface. A review of Schottky's derivation is found in Hutchison 
(Reference 6, pp 305 to 309). Although this section neglects space-charge 
effects, as did Schottky, space-charge effects will be brought into this 
theory and accounted for in Section IV. 

1. DERIVATION OF THE SCHOTTKY EQUATION 

As discussed by Ramo, Whinnery, and Van Duzier (Reference 22), the 
force exerted on an electron of charge-e located at position x away from a 
conducting plane is the same as the force exerted on the electron by a 
charge +e located at position -x. Then the electric field experienced by 
the electron is 

F = e/4ir< o (2x) 2 (129) 

The work done on an electron in bringing it from a great distance away 
from the conducting plane to position x is then 

x p 

W = -/ [(+e)(-e)/T6 ttf (x') ]dx 1 = -e(e/167rc 0 x) (130) 

CO 0 w 

If this work is defined in terms of a motive potential <f> so that 

m 

W = (-e) (131) 

then 

= e/16TiE o x (132) 
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This potential is shown in Figure 16 where the Fermi level is displaced 
above the vacuum level by an amount ■{■, the work function. Due to the 
way which <j> is defined, the reference zero potential is the vacuum 
level with no applied voltage. In principle. Equation 132 predicts that 
the motive potential is unbounded at the cathode surface so that the 
work function would have to be infinite. However, the idea of an image 
charge due to a perfectly conducting surface breaks down at distances of 
the order of interatomic spacings. As discussed by Nottingham, 

(Reference 13), this changes the shape of the motive potential very 
near the surface and leads to a finite work function. Beyond a few 
angstroms from the surface, the fact that the motive potential is 
accurately described by Equation 132 has been proven for many types of 
cathodes by the agreement of experimental data with the Schottky equation 
derived below (Reference 5, p 109). 


If a potential V is applied to a plane anode at a distance D from a 
plane cathode, the applied potential depresses the vacuum potential V x 
linearly with distance from the cathode surface according to 

V x = Fx (133) 


where 


F = V/D 


(134) 


is a constant electric field. The total field V{x), which is the sum of 
the fields from Equations 132 and 133, is expressed 


V{x> = e/16*r c Q x + Fx (135) 

and is plotted in Figure 17. The value of x for which V is minimum is 
found from the zero of the derivative to be 

*„,n * (,36 > 

The value of V at x„. is 
mi n 

<>37) 
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As seen in Figure 17, this represents a decrease in the work function 
from the value with no applied external voltage. The reduced work function 
is then 

$ r = <P - ^eF / /4nc Q (138) 

Using this modified work function in the Richardson equation. Equation 8, 
gives 

J = AT 2 exp {-U-/e'F/4ue )e/kT} (139) 

or 

J = J Q exp {(e/kT)/eF/4ne o } (140) 

where F = V/D. Calculating this effect for 5000 volts applied across 1 

millimeter gives x-„ = 84 A and V . = .085 volts or J = 2.3 J at 900°C. 

3 min min o Q 

Calculations for 1000 volts across 1 millimeter give x min = 190 A, 

V nn'n = ,038 volts and 0 = 1.4 J Q at 900°C. 

2. EDGE EFFECTS AND SCHOTTKY THEORY 

One goal of cathode testing is to determine basic physical properties 
of a cathode from specific electrical measurements. This allows cathode 
performance to be predicted independently of the test vehicle. For this 
purpose Langmuir space charge theory (References 7 and 8) and Schottky 
theory (Reference 4) are available to interpret electrical test data from 
the cathode. However, these theories are well developed only for certain 
simple geometries such as parallel planes and concentric cylinders and 
spheres. In practice the test vehicle only approximates one of these 
ideal geometries. When the theory is used tc interpret the measured data, 
errors are introduced. The calculated "basic properties" of the cathode 
are then dependent to a certain extent on the geometry of the test vehicle. 

A close-spaced diode which is often used for testing cathode emission 
is shown in Figure 18. Langmuir and Schottky theories for infinite planes 
are used to interpret the measured data. The most obvious error comes 
from electric field variations across the emitting surface due to electric 
field enhancement by the edge of the sleeve. Additional error comes from 


41 





AFWAL-TR-81-1019 


the finite extent of the emitting region which generally is smaller than 
the total cathode surface (see Figure 18). The presence of a boundary 
to the emitting layer is contrary to the conditions assumed in the 
derivation of Langmuir space charge theory for parallel planes. The 
space charge is dissipated more easily at the boundary resulting in 
greater current density than is predicted for the parallel plane case. 

This problem is discussed by Pierce (Reference 23) who shows that uniform 

emission can be obtained if a focus electrode is added at cathode 

potential which makes a 67.5 degree angle with the normal to the 

cathode surface (see also Reference 22). With or without the focus 

electrode the current varies with the three halves power of voltage under 
space charge limited conditions, as discussed by Langmuir (Reference 8, 
p 251) and Spangenberg (Reference 24). Without the focus electrode, the 
saturation is not abrupt due to the nonuniform current density. Also, 
the slope of the three halves power curve is not simply related to the 
square of the cathode-anode spacing. 

Focus electrode geometries other than the 67.5 degree sleeve can 
improve the uniformity of the space charge limited emission current. 
Computer programs for traveling wave tube gun design (Reference 25) are 
available which can evaluate more complex geometries. Although these 
designs lead to a more accurate representation of the emission by 
Langmuir space-charge theory for parallel planes, they will not improve 
the representation by Schottky theory under temperature limited 
conditions. 

Schottky theory for parallel planes becomes more applicable when an 
electrode at cathode potential (guard ring) extends the plane of the 
emitting surface. The guard ring gives a more uniform electric field 
once the applied voltage is high enough to decrease the space charge 
effects. Traveling wave tube gun design programs can also calculate the 
electric field for the test diode in the absence of space charge. 

Comparing the Schottky theory for the true field distribution with the 
Schottky theory for parallel planes shows the errors in the parallel plane 
assumption. 
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If the errors in the parallel plane assumption are not large, a 
simpler analytical approach gives an accurate approximation to the true 
field distribution. This approach begins with the calculation of the 
electric field about a plane edge above an infinite plane as shown in 
Figure 19. The Schwartz transformation relates the distance from the 
cathode edge to the electric potential and flux through a procedure 
outlined by Ramo, Whinnery and Van Duzier (Reference 22). 

The Schwartz transformation allows one to find a mapping from the 
complex z' plane into the complex z plane which preserves both the 
properties of the static electric field and the boundary conditions. If 
the z' plane is chosen as shown in Figure 20(a), the solution to the 
static electric field problem can be found in a straightforward manner. 
Let the boundary condition be V equals V Q to the left of x = 0 and V = 0 
to the right of x = 0. Then electric potential v and flux function u 
are given by 

u + iv = (V Q /tt ) 1 n {z 1 1 (141) 

where 

u = flux function =-/Edy+/Edx (142) 

x y 

i = v'TT 

F = electric field intensity 

A quick check of Equation 141 shows that 

v = (V Q /„) arctan {y’/x'l (143) 


satisfies both Lapace's equation and the boundary condition specified 
above as required for the electric potential. The inverse of Equation 141 


where 


z' = exp { ttw/ V q ) 
w = u + iv 


(144) 

(145) 


will be required below in the derivation. 
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For this particular problem the z plane would be that of Figure 19 
and the required transformation is (Reference 22, p 189) 

dz/dz' = K (z'-x' 1 /“l /T '" 1 ^(z'-x , 2 ) (ct 2 /lT ' 1) (146) 

The angles and u 2 and the positions x^ and x£ are shown in figure 20(a) 
and (b). Including the appropriate values for xj, x 2 , and in 
Equation 146 gives 

dz/dz' = (K/z’) /z^T (147) 

Integration of Equation 147 gives 

z = 2K[/z 7- T - arctan {/PCT}] + C (148) 

where C is the integration constant. The constant C can be evaluated by 
the requirement that the point z* - 1 transform into the point z = ih. 

This requirement imposed on Equation 148 gives 

C - ih (149) 

where h is the diode spacing. The constant K is evaluated from Eqjation 
147 through an integration from point B to point C shown in Figure 20 in 
both the z and z' planes. It must be that 

/ B C dz = / B C (dz/dz 1 ) dz’ (150) 

If the integration is about a small semicircle in the z' plane at B - C, 
then dz/dz 1 becomes arbitrarily close to iK/z’ and dz 1 can be expressed 
as iz'do 1 . Then Equation 150 can be expressed 

CSS dz = (ik/zoiz'de- (151) 

Evaluation of Equation 151 then gives 

K = i h/TT (152) 
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Finally, a combination of Equations 152, 149, 148, and 144 gives 

x + iy = (2ih/7i)[n/2 - arctan{/exp{7iw/V }’ - 1 } + /exp{ ttw/ V Q ) -1] 

(153) 

Since the arctangent is related to the logarithm by (Reference 19, p 80) 
arctan x = (-i/2) In {(1 + ix)/(1 - ix)> (154) 

the real part of Equation 153 can be expressed along the cathode 
surface as 

(x/h) = (1 /7i) 1 n{ (1 + /) - exp{ ttu/ V Q >)/ (1 - ST~^ exp{7ru/V o )) } (155) 

- (2/tt) ST - exp{7Tu/\y 

since w = u + io along the cathode surface. Note that u ranges from 
0 to - «> along the cathode surface. From Equation 142 the perpendicular 
component of the electric field F y is just du/dx. Taking this 
derivative in Equation 155 gives 

(h/V Q ) F y = d(u/V Q )/d(x/h) = -1/vT^Txyr^T} (156) 

The behavior of Equation 155 can be checked for u having large 
negative values. In this case Equation 155 can be approximated as 

x/h = (1 /ti) In {4 exp{-TTu/V Q }} - 2/tt (157) 

which says that 

u = -(V o /h)x (158) 

for large u and x. This is expected since u is just the line integral 
of the electric field. 
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The expression of Equation 155 is in an inconvenient form for 
further use since distance is given as a function of voltage. Attempts 
to get an analytical inverse have been unsuccessful. However, a 
numerical inverse is straightforward by means of the Newton iteration 
(Reference 20) also used in Subsection 5(d) of Section II. The Newton 
iteration is implemented for this problem by the function f{u} for a 
specified x such that 

fful = (1/tt) ln{ (1 + /I - exp{Ttu/V o lT/0 - A - exp{uu/V 0 ') J1 


— 2/tt) v 7 l - expUu/V '}' - x/h 

Then the desired u{x} is found by solving for 

flu) = 0 

by constructing the sequence 

U 1 = V /h + 1 ^ */ h > 

or 

U 1 = V o x/h if x/h < -2 

and 

u n + 1 = u n ' f{u n } /f '' u n } n = 1 * 2 » 3 


(159) 


(160) 


(161) 


(162) 


The choice of starting value u-j comes from an inspection of Equation 157. 
The function f{u p } = df/du comes by direct calculation on Equation 159 or 
by inspection of Equation 156 and is expressed by 

f {u n > = -/I - exp{TtiT/V q) (163) 


The series is truncated when hf{u}/x is sufficiently close to zero and 
this u n is taken as ulxl. Then F y {xl is calculated from Equation 156. 

If space charge is neglected and the surface electric field is known, 
the emission current density can be calculated from the Schottky 
Equation 140. A convenient procedure is to divide the cathode surface 
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into M rings of equal area, 
ring into equal areas is 

r- N 

The distance from the cathode edge to this mean radius is 

x^Nl = S + (d/2) (1 - /N/M - .5/M) (165) 

where 

S = sleeve width 

d = diameter of the emitting area 

A detailed sketch of this diode geometry is shown in Figure 18. The 

electric field F at a distance xIN; from the cathode edge is calculated 
y 

by means of the Newton iteration with Equations 155 and 156. This is, of 
course, an approximation since the cylindrical nature of the cathode is 
ignored in the field solution. The current density at each x{Nl is then 
calculated from Equation 140 and the average is taken over the M rings. 

The process is repeated for different values of applied voltage generating 
an I - V curve. The Newton iteration and Equations 155 and 156 have to be 
used only once to solve the field problem since the electric field is 
found in terms of the applied voltage. 

The solution for the current density and electric field for a 
particular example are shown in Figure 21. The parameters chosen are 
cathode diameter d = 3 mm, sleeve thickness S = 0.175 mm, cathode-anode 
spacing h - 1.5 mm, cathode temperature T = 900°C, and cathode-anode 
voltage V q = 1800 V. The voltage enhancement is 394 at the edge of the 
emission layer while the current enhancement is only about 8%. This is 
surprising at first glance, since the current varies as an exponential 
of the square root of the voltage. The coefficient in the exponent is 
small, however, so that the dominant function is the square root rather 
than the exponential for these experimental conditions. 


Then the mean radius which divides the Nth 
= (d/2)*'N/M - .5/M (1 


47 












AFWAL-TR-81-1019 


The average current as a function of voltage in terms of Schottky 
coordinates is shown in Figure 22. The effect of the edge is to 
increase slightly the average current density at a given applied voltage. 
Although the effect appears small in Figure 22, the increase in slope 
of the 1n{J/J Q } versus V curve gives an apparent decrease in temperature 
from 900°C to 834°C. The temperature is calculated from the Schottky 
plot by means of Equation 140 in the form 

T = /e/4 *.~ (e/k) ST /In{J/J s > (166) 

In this case the Schottky curve including edge effects is very nearly a 
straight line passing through J = J q at V = 0. The curve including 
edge effects is straight since the argument of the exponential in 
Equation 140 is small for this example. This gives support to the common 
practice of using the Schottky line intercept to find the zero field 
current even when the slope of the Schottky line doesn't agree with the 
measured temperature. 
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SECTION IV 

EXTENDING SPACE-CHARGE THEORY INTO THE ACCELERATING FIELD RANGE 

Figure I shows schematically the assumption made in the derivation 
of Equation 12. For this derivation, the potential due to the surface 
forces which the electrons must overcome occurs as a step function right 
at the cathode surface and has a value of ; electron volts below the 
Fermi level. Actually, the potential is not a true step function but 
acts over a short distance outside the surface due to the image force 
on the electron. Becker (Reference 5) points out that neglecting the 
finite extent of the image force has little effect on the space-charge 
distribution as long as there is a space-charge minimum volts below <t>. 
However, if there is no space charge minimum, the space charge potential 
interacts with the image force to lower the value of 4». 


1. BACKGROUND 

This interaction of the image force with the space charge occurs in 
much the same manner as discussed in the derivation of the Schottky effect 
in Subsection 1 of Section III. The derivation of the Schottky effect is 
based on the superposition of the image potential of the electron and a 
linear applied potential. The following discussion shows how the initial 
slope of the space-charge potential serves as the linear applied potential 
in the Schottky derivation. 

Due to its own image, an electron experiences a surface potential 

field 


V. = e/4'i (4x) 
image o 


067) 


where x is the distance from the cathode surface, and V. is measured 

i md 9“ 

relative to the potential a long distance from the cathode (vacuum level). 
In addition an electron experiences a space-charge potential field given by 
Equation 27 in normalized form. Figure 23(a) shows the image potential 
as a function of distance from the cathode surface while Figure 23(b) 
shows the space-charge potential as a function of distance from an 
example cathode. The example cathode of curve © has a current density 
J o of 5 A/cm 2 at 900 n C for dV/dx = 0 at the cathode surface. To illustrate 
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the behavior of the space charge potential for accelerating field 

conditions. Figure 23(b) also shows the potential-distance curve 

according to Equation 27 for a cathode with an accelerating field 
-4 

dV/dx = 6.44 x 10 volts/meter at the cathode surface. Although this 

2 

value of field strength (corresponding to the arbitrarily chosen 5.5 A/cm ) 

is chosen at this point for demonstration purposes only, the following 

discussion shows how this increased field strength gives an increased 

2 

current density of 5.5 A/cm at 900°C for the same example cathode. 

Note that the x scale in Figure 23(a) is expanded by a factor of 100 over 
the x scale in Figure 23(b). 

Figure 23 shows that over the range where the image potential is 

- 8 

effective (from the cathode surface out to about 5 x 10 meters) the 
space-charge potential can be approximated very accurately as an 
extension of the initial slope of the space charge curve 


V 

space 

charge 


{xl 


dV/dx 


(x)= F x 

cathode 

surface 


(168) 


Based on this assumption, namely, that the space-charge potential is 

linear near the cathode surface, the derivation for the Schottky effect 

applies if this F^ instead of V/D is used. The result is a potential 

minimum at a distance x . as in Equation 136 

min M 

x . = (e/F ) 1/2 /(16*r ) 1/2 (169) 

min ' x' ' o ' 


from the cathode surface, and the magnitude of the potential minimum 
is reduced from the work function 4> by an amount A<j> where from 
Equation 137 


A0 = (eF x ) 1/2 /(4nc 0 ) 1/2 


(170) 


For these conditions replacement of *j> with <|>-A<j> in Equation 8 gives 


J = AT 2 exp { -eU-(eF x /4TTc 0 ) 1/2 )/kT} 


(171) 
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Equation 171 with = 0 gives the zero-field current density 

J q = AT 2 exp { -e<>/kT} (172) 

which is the Richardson-Dushman equation. With this expression for J Q , 
Equation 171 becomes 

J = J Q exp{e(eF x /4TH_ o ) 1/2 /kT} (173) 


2. GENERATION OF J'VI IN THE ACCELERATING FIELD RANGE 

The previous discussion now provides a framework for calculating 
the voltage needed to give a specified current. Let V g be the voltage 
of the anode relative to the cathode. Then, referring to Figure 8, 


V 

app 


Vixl 


- V s + 


9 , - 


anode 


(174) 


where x = D - x c at the anode, D is the cathode-anode spacing, is 

the distance from the cathode to the potential minimum, and <j> and 6 are 

the work functions of the anode and cathode respectively. For simplicity, 

the remainder of this section takes = p and refers to the applied 

a c 

voltage as V. If J J , Equations 27 and 28 give the voltage according 
to the procedure outlined by Langmuir with n = 0. Equation 91 summarizes 
this procedure. 


If J > J , Equations 27 and 20 give the voltage 

V - (kT/e) n a (O/T /x 2 } (175) 

where n is calculated by a mixed Runge-Kutta-Taylor method using 

a 

n Q = (e/kT)(x 2 /v(J) F x (176) 

In Equation 27. The expression of fi 0 in Equation 176 is just an 
expression of the electric field in normalized units using Equations 87, 21, 
and 20. From Equations 173 and 176 

n Q = (4 1 r, o 3/2 /e)(kT/e) 7/4 (e/2 1 rm) 1/4 ln 2 fJ/J Q }//J (177) 

or 

h 0 = 6.4976 x 10' 2 T (7/4) ln 2 iJ/J q } //J 
with T in degrees Kelvin and J in amperes per square meter. 


( 178 ) 
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Figure 24 shows an example current-voltage plot for J ranging from 

2 

small values up through J>J Q . For this example J Q is 5 A/cm , T is 900°C 

and the cathode-anode spacing D is 1 mm. Figure 24 also shows the 
Child's law line. Equation 128, and the Schottky line. Equation 140 for 
comparison with this extended theory. For this example the extended theory 
matches Child's law closely up to J = 1.15 J Q but the approach to the 
Schottky line is only very gradual. 

3. BEHAVIOR OF J{V1 NEAR J EQUAL TO AND GREATER THAN J 

o 

At the point where the space-charge range meets the accelerating- 
field range, the space-charge minimum occurs right at the cathode. The 
derivative of Equation 175 with respect to J shows the behavior of the 
current-voltage curve from the accelerating-field side at the transition 
from the space-charge range to the accelerating-field range. Since n 
is also a function of J, the derivative is 


dV/dJ = (kT/e)[(D/2x, /J )d n /df, 

l . a 


f ' 

a 


+ (d n /dJ)d n /d h 

O d U 


where 


f ‘ = /J 0/x 0 

a l 


] (179) 


(180) 


Inspection of Equation 177 shows that dn Q /dJ = 0 at J = J Q . In addition 
an expression of n as an integral evaluated at r' along with the 
definition of the derivative shows that do /do is zero for o =0 

3 0 0 

(i.e., at J = J ). This expression for do,/d5 is 

0 3 0 

. i 

do /d5 = f [n ' 2 + 0 - P{*^1) explol - 0 - 24A)] 1/2 dC (181) 

a o ao Q 0 0 

The derivative evaluated at o Q = 0 is then 


do / dr'i 

a o 


o=0 


q[ + (1 - exp In) - (1 - 2v / o/iT)] 1/,2 dC 

'o o 


- / [(1 - Pf/^1) explol- (1 - 2/^)] 1/2 dc]/o 


(182) 
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For the left hand integral in Equation 182 choose same £ = A so that 
riff} satisfies 

\ = (1/2)[(1 - Pi/Ti) exp ini - (1 - 2/^h)] (183) 

Then a portion of the left hand integral can be bounded by 

/ 0 [i 1 q + (1 - PWn") )exp{n 1 - (1 - 2/n/V ) ^ df, < J Q [Hq]^ ^^ d£ = h Q A (184) 

Observe that the expression /l + x is bounded by 1 + x/2 for positive x 
since 

/T+ - * < .'1 + x + (x 2 /4) = 1 + x/2 (185) 


Use of this expression in the remaining portion of the left hand integral 
in Equation 182 gives 

r 1 f 1 

f [-ig+d-Pi vT7) )exp{r,}- (l-2 v / 777V) ] 1/2 dr. < / [(1 -Pi /T x ) )expf»,}-( 1 -2* nh) ] 1/2 d> 

A 0 

. \ 

f 

+ 1/2/ (r^/[(l-P{,n])exp|.,}-(l-2v^A)] 1/2 )d£ (186) 

o 

where the integrals on the right hand side have been extended to £ = 0. 
This preserves the inequality since the integrals are positive. Then 
from Equations 182, 184, and 186 the derivative is bounded by 


dn 


a 



lim 1 r- a + 

n ~ In A + 

n 0 n o 
o o 


r . 2 r,' 

-j- / [(l-P{/^l)exp{r l }-(l-2v / v7)]" 1/2 d£] 
0 


V° 


(187) 


A comparison of Figure 4 and Equation 183 shows that A must go to zero as 
n goes to zero. Then, since the dependence of n on h has to be 

0 a 0 

greater than or equal to zero. Equation 187 shows that 


dn,/dn 
a o 


= 0 


V° 


(188) 
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A combination of observations on Equations 177, 179, and 188 shows 
that the derivative of voltage with current at the zero field current 

J o 1s 


dV/dJ 


J=J o 


(kT/e)(D/2x 2 /j)dn a /df 


r = A 57 D/x 0 

o 2 


(189) 


which has a finite, non-zero, value. Since the slope at J = J Q is zero 
from the space-charge side and it has a finite, non-zero value from the 
accelerating field side, this theory predicts a discontinuity in the 
slope of J as a function of V at J = J . 

Figure 25 shows an example where this discontinuity in the slope is 
apparent. The dashed lines intersecting at J = J Q in Figure 25 show 
the limiting slopes as J approaches J q from the space-charge range and 
from the accelerating-field range. Van der Ziel (Reference 9) treated 
this example in the space-charge range, and his Figure 1 shows J 
approaching J y with a zero slope as does this Figure 25. Equation 189 
gives the limiting slope as J approaches from the accelerating-field 
range. Figure 25 shows that the limiting slope from the accelerating 
field range is almost the same as the slope of the Child's law line at 
J = J , although Child's law does not approximate J as a function of V 
for this example. The following discussion shows that the agreement 
between the limiting slope and the Child's law slope is not a coincidence. 

In many practical examples such as the one treated in Figure 24, 
the discontinuity in the slope of J{V} causes only a small perturbation 
in the curve. In this case the extended theory provides a nearly smooth 
transition from the space-charge range into the accelerating field range 
and then produces a gradual approach to the Schottky line. In this case 
the term dt c /dn in Equation 92 is small except for a very small region 
near = 0. Then Equation 92 reduces to Equation 189 and remains 
accurate for J somewhat greater than J . For many experimental conditions 
including those of Figure 24, Equation 27 which gives dn./dc , is 
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accurately represented by only one term - the one involving 2/nTA • 

a 

Consequently, Equation 189 reduces to 

dV/dJ = (kT/e)(D/2x 2 /J) 2 (n a M 1/4 (190) 

and n, in Equation 190 with their equivalent expressions 

a 

and 87 removes the normalizing factors and gives 
dV/dJ = (m/2e) 1/4 t 0 ' 1/2 DV 1/4 /J 1/2 (191) 

Rearrangement and integration of Equation 191 gives 

J = 4/9 , /2e/m V 3/2 /D 2 (192) 

which is the same as Equation 128, Child's law. This exercise shows 

that Child’s law is a good approximation for the current as a function of 

voltage not only in the space-charge range, but also a good approximation 

for a certain interval into the accelerating-field range. Figure 24 

hows, for that example, that J follows Child's law closely up to about 

J = 1.15 J . 

o 

4. FIELDS AND CHARGE IN THE INTERELECTRODE SPACE 

With the theory developed in Subsection 1 of Section II and Subsection 1 
of Section IV, it is now possible to plot the voltage V, the electric 
field F and the number-density of electrons N as a function of distance 
from the cathode for various conditions of cathode loading J/J Q . At a 
specified temperature T and current density J Q , a family of these plots 
can be made with d/J Q as a parameter. These curves show the changing 
conditions in the interelectrode space as the current drawn from the 
cathode changes, and they help to develop a firmer concept of the fields 
and charge in the interelectrode space. 

For space-charge-limited condition' Equations 28 gives the distance as 
a function of voltage to the cathode side of the potential minimum. After 
the addition of the normalizing constants from Equations 20, 21, and 87 
this relationship can be expressed as 

x = (x 2 //J) f, c {eV/kTl (193) 


Replacement of x 2 
from Equations 20 
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From Equations 84 and 85, V ranges from 0 to V $ in Equation 193 causing 

x to vary from 0 to x m with 
J m 

x m = f ' C tln{ V J} } (194) 

and 

V $ = (kT/e)ln{J o /Ji (195) 

The terms x and V are measured relative to the potential minimum; 
although it is more desirable to have x' and V' measured relative to the 
cathode surface where 

x' = x - x m (196) 

and 

V' = V - V s (197) 

Note that x is negative and V is positive while x' is positive and V' 
is negative in the cathode space. Then Equation 193 can be expressed as 

x' = (x 2 //j) c c < e (V' + V s )/kT} - x m (198) 

The electric field F is then -dV/dx where from Equation 28 

dV/dx = -(kT/e)(/J/x 9 )([1+P{/rO] exp{n 1 

1 C C (199) 

-0 * 3(n c /,) 1/2 ]) 1/2 

with 

n c = eV/kT (200) 

The number density of electrons from Equation 19 is 

N = ( E() /e) d 2 V/dx 2 (201) 

and from Equation 23 this is 

N = ( tQ kT/2e 2 )(J/x 2 2 ) exp{n c )[l + P{/^1] (202) 
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To the anode side of the potential minimum Equation 27 gives the 
voltage as a function of distance. After addition of the normalizing 
constants this relationship is expressed as 

V = (kT/e)n (x«tf /x 2 > 
or 

V' = (kT/e )n f(x‘ + x )/T /x 9 l -V 
a m 2 s 

The electric field is then -dV/dx where 

dV/dx = (kT/e) (v'J /xj([i - PirTT f]exp{ n J-[l - 2/rA) 1/2 

c. da a 

The number density of electrons from Equations 23 and 201 is then 

N = (, 0 kT/2e 2 )(J/x 2 2 ) expfn a )[l - Pf/^ 3] (206) 

For accelerating field conditions Equation 27 gives the voltage as 
a function of distance with the added condition on h Q from Equation 177. 
The expression for the voltage is then 

V = (kT/e)ii {x/J /x 9 ! (207) 

a C 

wi th 

n Q = (4 - 7e)> o 3/2 (kT/e) 7/4 (e/2T,m) 1/4 ln 2 {J/J Q l//j (208) 

Then from Equation 27 the electric field F is given by -dV/dx where 

1 /2 

dV/dx = (kT/e)(/J /x„)(n 2 + [l - Pl/O] exp{nj-[l - 2(n a M V2 ]) (209) 

C U a a a 

In addition, from Equations 23 and 201, the number density of electrons 
is given by 

N = ( !o kT/2e 2 )(J/x 2 2 ) explain - P(/^ 1] 
with the added condition on h in Equation 208. 


(203) 

(204) 

(205) 
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Figures 26 through 28 show the fields and charge in the interelectrode 

2 

space for an example cathode at 900°C with a current density of 5A/cm . 
Figure 26 shows the electric potential as a function of distance away 
from the cathode surface. The curve at J/J Q = .6 shows clearly the 
potential minimum. The curve at J/J Q = 1 shows the zero slope of V{x} 
at x = 0. The curves for J/J Q >1 show the linear behavior of V{x} near 
x = 0. The function J{V} can, in principle, be constructed from 
Figure 26 by drawing a line of constant x at the distance x = D represent¬ 
ing the diode spacing. The intersection of this line with the J{V} lines 
then gives the voltage associated with each value of J for the specified 
diode spacing. This line of constant x must fall to the right of the 
potential minimum to generate a valid J{V} curve, however, since the 
retarding field range requires a non-zero h 0 in the solution for Tin }. 

In contrast to the accelerating field field range, this A for the 
retarding field range depends on the spacing D and cannot be incorporated 
into the series of curves in Figure 26 which are independent of D 
(see Reference 9). 

Figure 27 shows the derivative dV/dx which is the negative of the 
electric field F. The curve at J/J = .6 shows -F to be negative near 
the cathode. The electric field goes to zero at the potential minimum, 
and then becomes positive in the anode space. The curve for J/J Q = 1 
shows F - 0 at x - 0 as required for the zero field current density. 

As J increases above J Q the curve for -F rises rapidly near the cathode 
and then levels off at distances away from the cathode. This tendency 
becomes more pronounced for large J/J Q until the electric field is nearly 
constant in the interelectrode space. This constant electric field is, 
of course, expected as space-charge effects become small. 

Figure 28 shows the number density N of electrons as a function of 
distance from the cathode. For small J/J Q the number density of electrons 
is large near the cathode due to electrons returning from the potential 
minimum. However the number density falls off rapidly away from the 
cathode toward the potential minimum. For J/J q closer to 1, the number 
density of electrons is smaller near the cathode but it does not fall off 
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so rapidly going away from the cathode. From Equations 195, 200, and 202 
the electron number density at the cathode is 

Nix = 01 = (« o kT/2e 2 )(J o /x 2 Z )[l + PilntJ Q /J>}] (211) 


At the zero field current density J , Equation 211 shows that 


Nix = 01| = ( i:o kT/2e 2 )(J o /x 2 2 ) 

J=J 


( 212 ) 


As J/Jq becomes very small. Equation 211 shows that the number density 
of electrons becomes twice that given in Equation 212, 


Nix = 01 


- (, 


0 kT/e 2 )(J 0 /x 2 2 ) = 


2N(x = 01 


J = 0 


(213) 


J=J. 


since the error function Plu! tends to 1 for large argument. Equation 213 
is reasonable since at J - J none of the electrons are returning while 
at J = 0 all the electrons are returning to give twice the density at the 
cathode surface. As J becomes larger than the electron density at the 
cathode is again larger than it is at J = and it falls away much more 

rapidly at distances away from the cathode. From Equation 210 the electron 
density at the cathode surface is 

N'x = 01 = (- Q kT/2e 2 )(J/x 2 2 ) (214) 

which increases linearly with increasing J. 

Figure 28 suggests that the electron density decreases away from the 
cathode with a finite slope for J less than or equal to J Q . However for 
J>J 0 the electron density has a decreasing slope with unbounded magnitude 
at the cathode surface. This can be shown rigorously from the equations 
involved. From Equations 28 and 202 

dN/dx = -(> 0 kT/2e 2 )(J 3/2 /x 2 3 )[exp{n c } (1 + Pl/^ }) (215) 

+ 1/ * / -n c ’][(l + P 11 )exp { n c } -(1 + 2/n^n )^ 
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1 


For large n this tends to the limiting value 

dN/dx = -(/2 c o kT/e 2 )(J o 3/2 /x 2 3 ) (216) 

which gives the slope at the cathode surface for J much less than J Q . 

Here use is made of Eguations 195 and 200 and it is noted that 

l//n c and /n^/expln } both become small for large n . It is also assumed 

for this illustration that the anode is sufficiently far away to allow a 

potential minimum in the cathode space. For n c near zero, corresponding 

to J near J , Equations 45, 54, and 215 show that dN/dx has the limiting 

value 

dN/dx = -(c Q kT/2/n e 2 ) (J Q 3/2 /x 2 3 ) (217) 


This can also be shown directly from Equation 215 by the use of L. 
Hospital's rule to show that 


1 im 


[1 + P{/n7l]exp{n }-[l + 2/n/Ti ] 

0 (I/2/77") ( -- c - - C - 2 - 


1/2 

) = (1/2/7) 


(218) 


In all these cases, then, from very small J up to J = J Q , the slope of 
the electron density dN/dx is finite at the cathode. 

For J>J q the derivative dN/dx is found from Equations 64 and 210 

dN/dx = (. o kT/2e 2 )(J 3/2 /x 2 3 )[(l - P{/^ }) exp ln a ) (219) 

-1 //tth" ][h ^ + (1 - P t/nTl) exp In) 

a u a a 

-1 ( 1 - 2/ri /71 ) ]^ 2 

a 

The term l/>nj increases without bound as n a approaches zero showing 
that the slope of the electron density is unbounded at the cathode. 

This is true for any fi Q >0 corresponding to J > d Q - 
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5. EFFECTS OF DIODE SPACING ON J■V( 

The zero-field current density J is often estimated from experi- 

o 

mental data through a determination of the "space-charge breakaway point." 
This is the point where the JiVl curve departs by a specified amount 
from Child's law. One method of specifying this "breakaway point" uses 
the projections of the slopes of the JiVs curve from the space-charge 
region and from the accelerating-field region plotted on two-thirds 
power paper. The intersection of the two lines representing these 
projected slopes is taken as an estimate of J . The dashed lines in 
Figure 29 give an example of this method. Affleck (Reference 17) 
presents experimental data showing that the diode spacing affects the 
value of the "breakaway point" determined in this manner. The following 
discussion uses the extended theory developed in this report to show how 
the "breakaway-point" technique leads to erroneous values of J q which are 
dependent on the diode spacing. 

The extended theory of Subsection 2 of Section IV is used to calculate 

the J(V• curve for a cathode similar to the one tested by Affleck. 

Figure 29 shows the current density as a function of applied voltage for 

two values of diode spacing, 0.3 mm and 0.9 mm. For the example the 

2 

temperature is 900 "C and the zero-field current density is 5 A/cm . The 

dashed lines in Figure 29 show that the "space-charge breakaway point" 

gives values of "J " higher than the true value of J in each case. The 
J o J 2 0 

projected slopes intersect at "J " - 7.6 A/crn for the 0.3 mm spacing and 
2 0 

at "J o " = 6.3 A/cm for the 0.9 mm spacing. This property of the 
"breakaway point" giving higher values of "J o " for the close-spaced diodes 
than for those with wider spacings agrees with Affleck’s conclusions 
based on data from a diode with variable electrode spacing. 

The above example shows that the "breakaway point" gives only a 
very crude estimate of the value for J . In the 0.3 mm case the estimate 
is 50% too high, while in the 0.9 mm case it is 25% too high. One reason 
for the erroneously high values of J Q is that J{V> follows the Child's law 
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line to higher values of J above J Q in the close-spaced diode than for 
the one with a wider spacing. A second reason is that the comparison 
for the cases with different diode spacing is made over the same range of 
diode voltage. Since the electric field is higher in the close-spaced 
diode the current density rises to a considerably higher values in this 
case. This causes the line projecting the slope of J1V) from the 
accelerating-field range to lie in the region of higher values of J for 
the close-spaced diode. In conclusion, the use of intercepting slopes 
on two-thirds power paper to determine the "space-charge breakaway point" 
leads to erroneous values of J Q which are dependent on the diode spacing. 
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SECTION V 
CONCLUSIONS 

A method for treating space-charge effects in the accelerating-field 
range is given. The method results in the current as a function of 
voltage for close-spaced, parallei-plane diodes. The space-charge theory 
developed by Langmuir is applied to the accelerating-field range along 
with the image-potential theory suggested by Schottky. The key to 
combining these two theories is to note that the space-charge potential 
is essentially a linear function of distance over the range of action of 
the image potential. The combined theory produces a current-voltage 
curve which begins at the zero-field current and voltage predicted by 
space-charge theory and then gradually approaches the Schottky line. 

For those conditions where Child's law approximates the current-voltage 
curve in the space-charge range. Child's law also approximates the 
current-voltage curve generated by the combined theory for an interval 
of current values extending above the zero-field current. 

The solution procedure outlined by Langmuir is summarized in a single 
equation. The derivative of this equation shows more explicitly the 
behavior of the current-voltage curve at the transition points from the 
exponential range to the space-charge range and from the space-charge 
range to the accelerating-field range. The current and the derivative of 
the current are both continuous functions of voltage at the transition 
from the exponential range into the space-charge range. At the transition 
from the space-charge range into the accelerating-field range the current 
is a continuous function of voltage, but the derivative of the current is 
discontinuous. This discontinuity arises because Langmuir space-charge 
theory always gives a zero slope for the current-voltage curve at the 
zero field current. However, the combined theory gives a finite, non-zero 
slope for the current-voltage curve at this transition point. 
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6 REM Data is dEta/dXi at X = 0 , delta X,print steps,Xmox 


15 

READ D0 , H , L , B 


20 

DSPLAY "DELTA Xi 

= " ,H,"dY^dX(at X = 0) 

22 

D S P L A Y 


30 

LET X = Y = Q = 0 


4 0 

LET E = 1 '' L 


90 

DSPLAY "Value of 

X i " , " V a Iu e of Eta" 

1 O0 

G 0 S UB 3 2 0 


120 

GOTO 160 


130 

LET U = Y 


1 3 1 

G 0 S U E 10 0 0 


1 ■-> /. 

LET F 1 = P 5 -(1-2 + 

U 1 s 1 .77245)+ D 0 ■# D 0 

136 

IF F 1 < 0 LET F1 = 

0 

137 

LET F1=SQR(F1) 


1 38 

LET U = Y + H/3 * F1 


139 

G 0 S U B 10 0O 


140 

LET Y=Y+H*(F1+H 

/ 4 * P 5 ) 

150 

LET X = X + H 


160 

IF (X+H-B)+1E-0 

7 >= 0 THEN 190 

170 

LET R = X + H 


130 

GOTO 200 


I 9 0 

LET R = B 


2O0 

LET fl=INTIE*R)/ 

E 

2 1 O 

IF ( A + L - R )-1E-0 

7 >= 0 THEN 230 

220 

LET A = A + L 


2 30 

IF ( X-A) +1E-07 

>= 0 THEN 260 

240 

LET G = A 


2 50 

G O S U B 3 2 0 


2 60 

IF R = E THEN 280 


270 

GOTO 130 


280 

IF ABS'R-AJ<.5* 

.000001 THEN 311 

290 

LET G = R 


300 

G 0S U E 3 2 O 


3 1 1 

STOP 


320 

LET U = Y 


32 1 

G 0 S U B 10 0 0 



322 LET F1=P5-(1-2*U1/1.77245)+D0»D0 

323 LET F1=SQR(F1) 

324 LET U = Y +( Q - X ) /3 * F 1 

325 LET X1= X 

326 LET X = Q 

327 GOSUE 1000 

3 2 8 LET Yl=Y+{0-Xn* (F1 + (Q-X1)/4*P5) 

329 LET X=X1 
370 DSPLAY 0 , Y1 

330 RETURN 

390 DATA 0,.01, .01, .2 

980 REM THIS SUBROUTINE CALCULATES (1-ERF(SQRCY))) *EXP (Y) 
1000 LET U 1 = S Q R ( U ) 

1010 LET T=l/fl+.327591*Ul) 

1020 LET P5 = T*« 1 .42141+T* (-1 ,45315+T*1 .06141 ) ) 

1030 LET P5=T*(.25483+T*(-.2845+P5)) 

1040 RETURN 
1040 RETURN 


Figure 6. BASIC Language Program Using Mixed Runge-Kutta-Taylor Method 
to Solve for n,U). 

d 
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0 EL T R xii 

Value of X i 
0 

.01 
. 02 
.03 
.04 
.05 
.06 
.07 
.03 
.09 
. 1 
. 1 1 
. 1 2 
. 1 3 
. 1 4 
. 1 5 
.16 
. 1 7 
. IS 
. 19 


.0f 

Value of Eta 
0 

.000025 
9.97337E-05 
2.23642E-04 
3.96303E-04 
6.19072E-04 
8.89S64E-04 
1 .20892E-03 
1 . 57603E-03 
1 . 9 9 1 1 6 E - 0 3 

2.4 5 3 7 6 E - 0 3 
2.9 6 3 6 6 E - 0 3 

3.5 2 G 6 2 E - 0 3 
4 . 12445E-03 
4.77433E-03 
5.47154E-03 
6.21435E-0 3 
.007003 

? .83726E-03 
8.716 9OE-03 
9.641 69E-03 


d V'7 d X 1 a t X = 01 


Figure 7. Sample Output from the Proaram of Figure 5. 
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Figure 8. Schematic of the Potential Distribution Between the Anode 
and Cathode under Space-Charge Limited Conditions. 
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J Q = 10 itiA/cm^ 
T = 700"C 
D = 0.2 nr.i 


d(J/J Q )/dV 


d 2 (J/J o )/dV 2 


APPLIED VOLTAGE V 


LVJLIs J 


Figure 9. The Current Density and its First and Second Derivatives as a 
Function of the Applied Voltage for an Example Cathode, 
for this example : is assumed equal to * . The dashed lines 
cover the retarding field range while the c continuous curves 
cover the spare charge range. The value of the voltage at 
.1 J is 0.9994 only by coincidence in this particular example 
















[A/m 4 ] 
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x = 4.58 
T = 900°C 
D = 1.0 mm 


x = 1.05 
T = 700°C 
D = 0.2 mm 


.3105 

LIMITING 

VALUE 


[A/m 4 -] 


5.864 
LIMITING 
VALUE 


Figure 11. The zero field current J Q as a Function of J g the Current at 
the Transition from Retarding Field to Space-Charge-Limited 
Operation for D = lmm, T = 900 C and for [) = 0.2 mm, T = 700°C. 









2 D S P L ftY "Program of 3 Nov 30" 

4 DS P L ft Y "Prints J < as a function of Jo" 

5 D S P L ft V 


10 

LET 

T = 

10 73 

20 

LET 

D = 

. 0 0 0 1 

30 

D S P L ft V 

" T E H P = ( ' K ) " , T , “ D I S T = ( M ) " , D 

45 

LET 

ft 2 

= CT 1.751*1.08 8 6E-06 

49 

LET 

N = 

. 1 

50 

D S P L ft V 

" D X 2 = " , D / ft 2 

5 1 

D S P L ft Y 


52 

D S P L ft Y 

" J o " , " J < " , " V < " , " X I c " 

5 3 

FOR 

0 = 

1 TO 4 STEP 1 

54 

LET 

N = 

N * 1 O 

5 5 

FOR 

I 5 

= N TO 10 * N STEP N 

60 

LET 

I = 

6.5 * i: ft 2 - D i * ( ft 2 /" B ) 

65 

I F 

I > 

= 15 LET I= .5*15 

70 

LET 

E 1 

= - L 0 G t I I 5 I 

75 

I F 

E 1 < 

O LET E1= O 

9 0 

LET 

V 1 

=El*T*8.62E-05 


1O0 LET E = E 1 

110 GoSU £ 14 0 0 

12 0 LET H = Z + D "ft 2 * $ QR( I 1 

2 3 0 LET T 1 = 1 i 1 + .327591 * S Q R (El 1 ) 

29 0 LET P 5 = 1 .42 141+(-1 .45315+1 ,06141*T1)*T1 
291 LET P 5 = ( .25 433+ i - .284497+P5+T1 ) * T 1 ) * T 1 

3 0 0 LET F 2 = 2 ■* E X P (El ) - P 5 - 1 - S Q R (El ) * 1 . 12338 

3 10 IP F 2 <= 0 LET F2=1E-12 

3 2 0 LET F 2 =1/3 0Rf F 2) 

322 LET H 1 = F2- I + Bs2 ' R2/SQR( I ) 

3 2 4 LET I 1 = I -H - H1 

3 2 6 IF ft E S ( t I - I 1 i r-1) < . 0 0 0 1 GOTO 332 

323 LET 1=11 
3 3 0 GOTO 7 0 

3 3 2 D S P L ft V I 5 , I , V 1 , Z 

350 NEXT 15 

351 NEXT 0 
1336 STOP 

1390 REM THIS SUBROUTINE CRLCULRTES ZC RS fl FUNCTION OF ETRC 
14 0 0 IF E > = .2 THEN 1450 

1410 LET RsSGR(E) 

1420 LET Z=-2*R+.376126*E+2.51956E-02*E*R+6.9163E-04*E*E 

1421 LET Z = Z- 9.9 4 O 9 E-0 4 * E * E * R 
1440 RETURN 

14 5 0 IF E > = 3 THEN 1490 

14 6 0 LET Z = -( - 6.15614+E-65.6912«E*E-59.4893*Et3+4.10335*Et4) 
1461 LET Z=Z/262909-14.6202*E-49.5951*Et2-16.0639*Et3+Ef4) 
1480 RETURN 

1490 LET Z=-2.55389+1.41421*EXP(-E/2)-.0123*EXP(-E) 

1491 LET ZsZ+1/3^1.41421*(SQR(E/3.14159)+l)*EXP(-3*E/2) 

,rio RC 7 

i6<ro 

Figure 12. BASIC Program to Solve Equation 113 and Generate J g and V g 
as Functions of J Q . 
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,e 1 . I he /pro field Current .1 and the Transition Current J c Plotted 

.1 . turn lions of the Transition Voltage V for 0 0.?mm, T - 700 C. 

I 


ZERO-FIELD CURRENT [A/ni' 








r.V e /kT [NORMALIZED UNITS] 

Figure 14. Graphs of J p and J p as Functions of V . These Universal Curves 
Relate the Measured Current at the Exponential Point to the 
Zero field Current of the Cathode. 
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external field is 10' volts/meter. 
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DIODE 

SPACING ANODE 

////i/ / // // // // // // / u_t /J 

—|— j v:-rjr_ - rn CATHODE 

i ! 

-' r .- SLEEVE INNER DIAMETER 

—• ! J — SLEEVE OUTER DIAMETER 


Figure 18. Schematic of a Close-Spaced Diode Tester Showing the 

Cross Section of a Cylindrical Cathode and Anode. The 
diameter of the emitting layer is the same as the sleeve 
inner diameter. 


U-0 

y 


* 


T 


v=o 


77777777777777777777777777777777 777 777 7 T7T7T71 v=v o 

x 


Figure 19. A Plane Edge Above a Conducting Plane. The flux function u 
is set to zero at X = 0, y = h. 
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x/h NORMALIZED DISTANCE 


Figure 21. Current Density and Electric Field at the Cathode Surface as 
a Function of Distance from the Cathode Edge - All in 
Normalized Units. Here S = sleeve thickness, h = cathode to 
anode spacing, V Q = applied voltage, T = temperature of the 

cathode, d = diameter of the emitting area, X = distance from 
cathode edge, and J = zero-field current density. 
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Figure 22. Theoretical Schottky Plots. The curve "with edge effects" 
takes into account the enhanced field at the cathode edge. 
The curve "without edge effects" refers to infinite parallel 
planes. 



V POTENTIAL (VOLTS) 












J/J CATHODE CURRENT (NORMALIZED) 
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Figure 25. Current-Voltage Plot Showing a Break in the Slope of J as a 
Function of V at J = J . The dashed lines show the limiting 

values of the slope at J = J . The Child's law line and the 
Schottky line are shown for comparison with the combined theory 
developed in this paper. 
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